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Abstract. Many recent studies have highlighted certain failures of the standard Eulerian- 
space cosmological perturbation theory (SPT). Its problems include (1) not capturing large- 
scale bulk flows [leading to an 0(1) error in the 1-loop SPT prediction for the baryon acoustic 
peak in the correlation function], (2) assuming that the Universe behaves as a pressureless, 
inviscid fluid, and (3) treating fluctuations on scales that are non-perturbative as if they 
were. Recent studies have highlighted the successes of perturbation theory in Lagrangian 
space or theories that solve equations for the effective dynamics of smoothed fields. Both 
approaches mitigate some or all of the aforementioned issues with SPT. We discuss these 
physical developments by specializing to the simplified ID case of gravitationally interacting 
sheets, which allows us to substantially reduces the analytic overhead and still (as we show) 
maintain many of the same behaviors as in 3D. In ID, linear-order Lagrangian perturbation 
theory (“the Zeldovich approximation”) is exact up to shell crossing, and we prove that n**^- 
order Eulerian perturbation theory converges to the Zeldovich approximation as n —)• oo. In 
no ID cosmology that we consider (including a CDM-like case and power-law models) do 
these theories describe accurately the matter power spectrum on any mildly nonlinear scale. 
We find that theories based on effective equations are much more successful at describing the 
dynamics. Finally, we discuss many topics that have recently appeared in the perturbation 
theory literature such as beat coupling, the shift and smearing of the baryon acoustic oscil¬ 
lation feature, and the advantages of Fourier versus configuration space. Our simplified ID 
case serves as an intuitive review of these perturbation theory results. 
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1 Introduction 

“[One] grows stale if he works all the time on insoluble problems, and a trip to the beautiful 
world of one dimension will refresh his imagination better than a dose of LSD. ” 

Ereeman Dyson 

Observations of the large-scale matter distribution show structure on a wide variety of 
scales. The organization of galaxies into a beaded, hlamentary “cosmic web” appears to be 
the natural outcome of gravitational instability in a cold dark matter dominated universe 
acting on an almost scale-invariant, Gaussian initial matter density field [1-3]. The process 
of structure formation is a competition between gravity and universal expansion and, as such, 
is a powerful probe of the cosmology [4, 5]. Owing to the spectrum of matter fluctuations in 
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the Universe, smaller scale matter fluctuations on average become nonlinear, collapse, and 
virialize at earlier times than larger structures, with nonlinear evolution happening coherently 
on < 10 Mpc scales at present. The processes that drive the evolution on these scales can 
only be studied in detail numerically. However, when “smoothing” the cosmic density and 
velocity field on scales of > 10 Mpc, perturbative methods may still provide an accurate 
description of structure formation. 

Linear-order Eulerian cosmological perturbation theory has been tremendously useful 
for understanding the anisotropies in the Cosmic Microwave Background [6, 7] and the statis¬ 
tics of present-day matter fluctuations smoothed over ;§> 10 Mpc scales [8]. The standard 
formulation of post-Recombination, Eulerian cosmological perturbation theory solves the 
continuity and Euler equations, treating the matter as a pressureless fluid, and solutions 
have been found to all orders in the initial matter overdensity [9-11]. Zeldovich [1] proposed 
a different perturbative approach that, instead of solving for densities and velocities, solves 
for the matter displacement field that is consistent with linear Eulerian perturbation theory, 
the “Zeldovich approximation”. These displacements can then be used to predict the density 
field, providing a better description than linear Eulerian theory into the mildly nonlinear 
regime [Ij. Zeldovich’s Lagrangian theory has been extended so that the displacement can 
be calculated beyond linear order [12-15]. Eor a review of standard perturbative methods 
in cosmology see [16]. In the decades since Eulerian perturbation theory and Lagrangian 
perturbation theory were devised, the concentration has predominantly been on schemes to 
accelerate the rather slow convergence of these original theories on mildly nonlinear scales 
by resuming parts of the expansion [17-20]. However, resummation schemes do not address 
all of the deficiencies of the original perturbation theories, particularly that their solutions 
are only valid prior to shell crossing and that they treat nonlinear scales as if they were 
perturbative. 

Recently there has been renewed effort to address these deficiencies by starting with 
effective equations that depend only on overdensity fields that have been smoothed so that 
their RMS is less than unity and, hence, perturbative methods can be rigorously applied [21- 
25]. Indeed, [22] claimed that these effective theories at 2-loop order remain accurate at the 
percent level to even k k, 0.6 /iMpc“^, or ~ 4^ times more modes than where Eulerian pertur¬ 
bation theory (at the low orders that are calculable) does. A second recent development has 
been the construction of higher-order Lagrangian theories [26-28], following the realization 
that Lagrangian perturbation theories fare better at describing the large-scale advection of 
matter. Such advection leads to virtually all of the nonlinear evolution of the baryon acoustic 
oscillations (BAO) seen in the matter correlation function [29, 30]. Neither of these develop¬ 
ments has come without controversy. Lagrangian theories still make the same uncontrolled 
approximations as standard Eulerian perturbation theory. Current calculations in effective 
theories rely on arguments that reduce the number of diagrams/coefficients [24, 31, 32]. The 
remaining coefficients are often treated as free parameters that are fit to the measured power. 
It is of debate whether the successes of the effective theories relative to previous theories in 
matching the nonlinear power spectrum owe to this extra freedom. (The counter-argument, 
which we ascribe to, is that these extra parameters are required for a consistent theory. How¬ 
ever, there is no guarantee that the low-order expansions can be used reliably close to the 
non-linear scale, the limit most applications of perturbative methods are most interested in 
describing.) It may even be the case that a fully perturbative approach to describe mildly 
nonlinear scales will never be successful. A fully perturbative approach appears to be in 
conflict with (popular) halo models for large-scale structure - which have had success at 


- 2 - 


describing the nonlinear evolntion in the matter field [33-36]. In halo models, the beyond- 
linear-order evolution is largely described by the density profile of halos, which arises from 
very nonlinear processes that presumably are not captured by any perturbation theory. This 
intuition has inspired the recent development of hybrid models that use both perturbation 
theory and assumptions about the matter distribution in and around halos [37, 38]. 

In this paper, we break from the standard approach of solving cosmological perturbation 
theory in the full three dimensions, and instead consider it in one spatial dimension (ID). 
One spatial dimension corresponds to modes all oriented in the same direction - breaking the 
statistical isotropy generally assumed. Equivalently, it corresponds to the interaction of infi¬ 
nite sheets of matter where the force is independent of distance from each sheet. Computing 
the total gravitational force thus amounts to counting the number of sheets to the left and 
right. In our cosmological context, these sheets are moving in a Hubble flow relative to one 
another (and the dimensions traverse to the sheets are also in the same Hubble flow such that 
the surface density in each sheet scales as a“^). Despite this high degree of symmetry, the 
same dynamical equations (such as the continuity and Euler equation in the standard Eule- 
rian formulation) apply as in 3D and, by extension, ID cosmological perturbation theories 
make all of the same assumptions.^ However, both the numerical simulations and analytic 
calculations are significantly simplified in ID. The three-dimensional d^k integrals that occur 
in all 3D calculations collapse to one-dimensional dk integrals, allowing the computation of 
higher-order solutions more easily. Simulations are able to have much higher dynamic range 
in wavenumber in ID than in 3D for the same memory and operation count. [We show in 
Appendix A that to simulate the matter power spectrum to 1% within a wavenumber range 
of Ak = 0.01 Mpc“^ requires solving for the dynamics of ~ 10^ sheets, whereas to do the 
same at a wavenumber of A: ~ 0.1 Mpc“^ - a mildly nonlinear scale at z = 0 in our Universe 
- with Ak = 0.01 Mpc“^ requires almost 10^*^ particles.] The reduced cost of simulations 
allows us to test most of the assumptions of different perturbative approaches on a wide 
range of cosmological models. A final (and very important) advantage is that ID allows us 
to calculate the results of both Eulerian and Lagrangian perturbation theory at infinite order 
(both yielding the Zeldovich approximation). Intriguingly, we show that standard Eulerian 
and Lagrangian perturbation theories evaluated at infinite order do not yield a correct predic¬ 
tion for the matter power spectrum at any mildly nonlinear scale in any cosmology that we 
consider because these theories err at describing the dynamics around collapsed structures. 

Of course the Universe has three spatial dimensions, so working in ID loses some crucial 
aspects of structure formation. However, while the gravitational dynamics in ID differ from 
those in 3D, ID still shares many features with 3D. Most trivially, in the limit that the 
wavenumbers of two modes have much different magnitudes, their coupling has to be the 
same (up to a geometric factor) in ID and 3D. Less trivially, we show that many results that 
have recently raised excitement in the perturbation theory literature, such as beat coupling 
[39], the smearing and shifting of the BAO peak owing respectively to the RMS matter 
displacement [28, 30] and the coupling to large-scale modes [40], and the (un)importance 
of stochastic terms in the nonlinear evolution [31, 37, 38], have analogues in the ID case. 
Indeed, we find that discussing all of these recent developments in our simplified ID case 
serves as a simple and intuitive review of these perturbation theory results. 

All of the calculations in this paper will assume Newtonian gravity, similar to most 
previous beyond-linear order perturbation theory studies. This approach is traditionally 

^The only assumption that ID does not make that underlies almost all 3D perturbation theories is that 
the nonlinear velocity field is curl-free. 
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Figure 1. Comparison of the correlation function in the concordance ACDM cosmology (top panel) 
with the correlation function in our ID CDM-like cosmology (bottom panel). Both the ID and 3D 
cosmologies share the same linear dimensionless matter power spectrum. All curves are for 2 ; = 
0 except the green dashed curve, which shows the Zeldovich approximation at 2 : = 1 in the ID 
cosmology. In both ID and 3D, the Zeldovich approximation fares excellently at describing the 
nonlinear development on > 10 Mpc scales. This success can be gleaned in the ID panel by comparing 
the Zeldovich Approximation curve with the non-linear curve, which was computed with our iV-body 
code. 


justified because our cosmology is in the limit where the nonlinear scale is much smaller 
than the horizon scale - the scale at which General Relativistic corrections apply. Indeed, 
the Newtonian calculation correspond to a gauge choice in General Relativity in this limit 
[41, 42]. As it will further simplify our analysis, we shall focus throughout on an Einstein de 
Sitter background cosmology (fl^ = 1), commenting on the transition to other cosmologies 
where appropriate. Our calculations in the main body of the text are in comoving coordinates, 
X. We use u to denote the peculiar velocity, and V a derivative with respect to x. Finally, 
most of our expressions will not explicitly show the temporal argument for densities, power 
spectra, etc. 

This paper is organized as follows. In section 2 we introduce our notation and con¬ 
ventions, then describe the standard Eulerian and Lagrangian formulations of cosmological 
perturbation theory (sections 2.1 and 2.2), proving their equivalence, and lastly describe the 
effective field theory of large-scale structure (section 2.4), working out examples for a ID 
GDM-like cosmology. We break from this cosmology in Section 3, extending our discussion 
to ID cosmologies with power-law initial matter power spectra. Finally, section 4 discusses 
the lowest order corrections to the matter correlation function and power spectrum in all of 
these theories; we show that our ID setting allows many results of perturbation theory to be 
more simply derived. Some of the technical details are relegated to a series of appendices. 
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2 Standard perturbation theories in ID 

We define the matter overdensity power spectrum, -P(k) as P(k) = x (27r)^/F, where V 

is the volume and S(x) = p(x)/p—l, p(x) is the matter density at configuration space position 
X, and tildes indicate Fourier space using the convention that is standard in cosmology in 
which the 27r’s only appear under the dk’s. We specialize to the case in which modes are 
directed along a single axis such that P(k) can be written as 

F(kii,kr,z) = Pid;.(A:||)(27r)2<5^(k^), (2.1) 

and we define = A:||P(A:||, k_i_; 2;)/7r to be the dimensionless power spectrum. Throughout 
we shall assume a ID version of isotropy, “reflection invariance”, so that Pid(—k) = Pid(k). 

We will use a few different parameterizations for the linear-theory power spectrum in 
order to explore how different perturbative theories fare. Mostly we will choose a ID linear 
perturbation theory power spectrum, PL{k-, z), to have a CDM-like form with 

7T-^kPL{k) = {27rY^ePcDM{k; z), (2.2) 

where Pcdm is the linear theory CDM power spectrum calculated with cosmological transfer 
function codes such as CAMB,^. Here, as in what follows, we have dropped the argument 
for the perpendicular wavevector and the redshift in Pl. This CDM-like form results in the 
same variance in the density per interval in k as in 3D CDM, as well as the same linear-order 
parallel RMS displacement. The top panel in figure 1 shows the linear correlation function 
in the 3D case, and the bottom panel shows this statistic in the ID case. Also shown is 
the Zeldovich approximation estimate for the nonlinear evolution (which works excellently 
in both cases; [28, 30]; section 2.2.1). In ID, the BAO feature is somewhat sharper, and also 
the nonlinear evolution at 2 ; = 0 is more substantial. Empirically, we find that the nonlinear 
evolution is more similar if we compare the ID solution at z = 1 to the 3D one at z = 0 
(compare the thin dashed curve in the bottom panel to the thick dashed in the top panel of 
figure 1), an identification used later on. 

In section 3, we further consider power-law linear-theory power spectra, parametrized 
as 

{2^)-^kPdk) = , (2.3) 

for A; > 0 and the symmetric form for A: < 0. The power-law case is particularly convenient 
because it admits self-similar solutions in our Einstein de Sitter cosmology, allowing the 
mapping of its nonlinear statistical state at one epoch to that at any other epoch. 

2.1 Standard (Eulerian) perturbation theory 

Standard [Eulerian] perturbation theory (SPT) makes the approximation that the Universe 
behaves as an ideal fluid, without pressure or viscosity, solving the continuity and pressureless 
Euler equations in the presence of gravity [9-1 1, 16, 43]. In ID these equations are given by 

d^8 + e = -V{6u), (2.4) 

drO + ^9 + 4TrGa^p6 = -V(uVu), (2.5) 

^http://camb.info 
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where the latter equation is the gradient of the Euler equation with 6 = Vn, dr = dt/a, 
% = aH^ and V denotes the derivative with respect to the comoving coordinate. Our ID 
symmetry has allowed us to make u a scalar. In Fourier space eqs. (2.4) and (2.5) become, 


~ ~ dk' k ~ ~ 

dr6{k) + e{k) = - —- 9{k') S{k - k'), ( 2 . 6 ) 

J—oo ^ 

~ ~ ~ r°° rik' P ^ ^ 

drOik) + n e{k) + 47rGa^p6{k) = -J ^ ^ 2k'ik-k') ~ 

These equations can be solved perturbatively with the ordering 

OO OO 

5(fc,t) = ^ a^{t)6^'^\k), e{k,t) = -n{a) ^ a^{t)6^^\k), (2.8) 

m=l m=l 


in an Einstein de-Sitter Universe (as considered here), where a is the scale factor and is 
the order solution for the overdensity. (It is common practice to replace o'” —>• D"*, where 
D is the linear theory growth factor, in order to treat cosmologies that are not Einstein de 
Sitter. This replacement is actually exact in ID.) In addition to assuming that the universe 
is a perfect fluid, the SPT expansion assumes that even though (5(x) > 1 will often apply, 
the perturbative solution to these equations is still valid when smoothing the solution on a 
scale R such that {d)R < 1. 

The solution at each order can be written as 


/ Sj^XkA 


= a 


dki....dkn D 

(27r)»"-i 


ki - k) 


( Fniih }) 

\ Gn {{ ki }) 


Sliki) ■ --SLikn), 


(2.9) 


where we have unconventionally kept the expression in terms of the linear overdensity at a, 
6l, rather than at o = 1, and the and Gn obey recurrence relations [10, 11, 16]: 


Fn — (2t1 + l)Xn + Yn ’, Gn — 3Xn + uYn - 


The argument of each function in these relations is (ki, - ■ ■ ,kn), and 


= 


Yn = 


1 ^ }Z 

77, I yy? yy / , 'r^Gmiki, ...^km)Fn—m(km+i, ■■■ikm), 
(Zn + 6)(n — 1) Ki 

m=l 

1 

To I TT / > TZ Gm{kl, ■■■km)Gn—m{km+lT ■■•T km), 

[Zn + 3)(n - 1) K 1 K 2 

^ m=l 


( 2 . 10 ) 


( 2 . 11 ) 

( 2 . 12 ) 


with Ki = ki + ... + km, K 2 = km+i + ■■■ + kn, k = Ki + K 2 , and Fi = Gi = 1. For 
example, the recurrence relations yield F 2 {k,k') = G 2 {k,k') = 1 + [k/k' + k'jk)j2 and 
F^'^[k,k',—k') = A:^/(6P^), where superscript ‘sym’ indicates the symmetrized form - Fn 
averaged over all permutations of its arguments - which is the form that is easiest to use in 
calculations. Using these solutions for the Y'^\k) to fourth order in 5^ (also referred to as 
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1-loop’) the matter power spectrum is given by 


pl-loop/, N 
-TgpT 

P22{k) 


Pi3{k) 


Pll + P 22 + -P 13 , 
/■“ dk' 


^Frik',k-kr 


dk' 


k-k' 


+ 


/ °° rlk' 

^Frik,k', 


— k‘^r]‘^PL{k), where = 


PL{k')PL{k-k'), 

^PL{k')PL{k-k'), 

-k')PL{k'), 

/■“ dk' PL{k') 
y_oo 27r A:'2 


(2.13) 

(2.14) 

(2.15) 

(2.16) 
(2.17) 


and we have adopted the notation Pnm{k){2TT)6^{k — k') = {k')) 

The 1-loop power spectrum in SPT, Pgp.|?°^(A;), receives some contribution from modes 
with k' 3> k, modes that in many cosmologies are highly nonlinear. That SPT shows any 
sensitivity to nonlinear modes motivates the effective theories discussed in section 2.4. In the 
opposite limit in which /c' <C A:, expanding Pi^{k — k') in /c', the leading dependence is^ 


4 pt » ^') = Pdk) - 2kVkPL{k) [ ^PL{k'). ( 2 . 20 ) 

4|fc'|<|fc| 2vr 


Thus, the sensitivity of Pspt(^) to a long wavelength mode {k' <C k) is via the coupling of 
its gradient to the large-scale variance, which acts to whiten the power spectrum and smear 
out features. The same coupling to large-scale modes holds in the Lagrangian and effective 
theories discussed later (as all theories discussed here have the same infrared behavior). 

The coupling of P{k) to a long wavelength modes, as given by eq. (2.20), is different 
from how the overdensity, 5, couples to such modes. The coupling of 5 to long-wavelengths 
modes, and the induced effects on the power spectrum, is often called ‘beat coupling’ [39, 46] 
or ‘super-sample covariance’ [47]. It has been studied in [39, 47, 48]. The coupling of the 
small-scale mode 6{k) with linear overdensity diik) to the overdensity on some much larger 
scale 5v is 5{k) = [1 -|- 25v]5L{k) in ID, which differs only slightly from the 3D result: 
5(k) = [1 -|- 34(5y/21]5L(k) [16, 46, 49, 50]. To derive this result note that there are two 
pieces of the k' integral where the 2“'^ order contribution to 5{k) involves long-wavelength 

^In deriving eq. (2.15) we have simplified the P 22 term using the fact that within the convolution integral 
we can interchange the arguments of F 2 . Beginning with 


Fi(ki,k2) 



k2 ^ ki 


= 1-t 





(2.18) 


and using the exchange symmetry 


1/p fey 

4 Vfe fey ^ 2 2 'vfey 


(2.19) 


in the integral we have —>■ (3/2) + 2{k\/k2) + (l/2)(fci/fe)^. For this way of writing the integral, all of 
the IR divergences occur for k' ~ 0, rather than having some also occur at fc ~ k', as occurs on most forms 
of SPT expressions. A similar organization is possible in 3D [44, 45]. 

^This is the same as the more often seen k^rf PL{k) form. The difference comes because we have written 
eq. (2.15) to have all infrared divergences occur for k' ~ 0, rather than having some occur at A: ~ k', which 
leads to k^rj^PL{k) being exactly cancelled. 
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modes, (5i(e) (either \k'\ ~ e or \k—k'\ ~ e) and by symmetry both give the same contribution. 
Assuming F 2 and 5l are smooth and can be taken out of the integral we thus have 

/ +<: Ah' _ 

— 5{k'). (2.21) 

Now noting that the fc'-integral is just Sy = I-l /2 dx5{x) for large L and [A 2 (—e, k) + 

F 2 {+€,k)] —)• 2 as e —)■ 0, we obtain the desired result. This coupling is important for 
understanding the shift of the BAO peak (section 4). 


2.2 Standard Lagrangian perturbation theory 

SPT perturbatively expands the RMS relative displacement between points separated by r, 
which at linear order and ID is oc J" dk P{k)/k‘^(l — cos[A:x]). It so happens that in our 3D 
Universe and also the ID CDM-like cosmology the RMS displacement is the largest pertur¬ 
bative effect at X > 10 Mpc in the matter correlation function. Indeed its size, ~ 15 Mpc at 
X ~ 100 Mpc and z = 0, is comparable to the width of the BAO peak. The displacement 
originates primarily from infrared scales at which linear theory is a good approximation, sug¬ 
gesting that it should be fully calculable and should not be treated perturbatively as in SPT. 
Lagrangian perturbation theory does not expand in the relative displacement and, hence, 
has been found to fare much better at capturing the BAO peak in the correlation function 
[20, 26-30, 51-53, and section 4]. 

In its standard form, Lagrangian perturbation theory (LPT) attempts to solve pertur¬ 
batively the equation [16]®: 

T(g) + 2R'^(g) =-V())(g + T). (2.22) 


Similar to SPT, LPT solves this equation in powers of the linear density field. Aside from not 
expanding in the linear theory displacement, LPT has the same problems as SPT. Like SPT, 
LPT makes the assumption that even though 5 1 on certain scales the solutions are still 

meaningful when smoothed over a scale R such that {d)ji <C 1. LPT is also not valid through 
shell crossing as infinities appear that invalidate the perturbative expansion as shown below. 

The solution to eq. (2.22) can be used to calculate the real and Fourier space density 
fields via 


1 + <5LPT(a;) = J dq 5^[x - q - 'l'(g')] = det[l -|- VqT] 


-1 


x=q+'i>{q) 


dk 


= dq I — e 


^ik[x—q—^{q)] 

2tt 


4pt(A:) = J dq - l) . 

It follows that the power spectrum of the density is 
PhPTik) = L~^ 

= f dq (^e^N^di^m{lkAp)c{q) _ ^ 


dqidq2 (^(|g-*U’l'(<?i)-^(g2)]^ _ i 


(2.23) 

(2.24) 

(2.25) 

(2.26) 
(2.27) 


® As with SPT, this approach in 3D assumes vorticity is zero and hence only considers the equation for the 
scalar part of the displacement, a technicality of no importance to our ID analysis. 
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where A = 'l'(gi) — '£'( 92 ) and L is the integration “volume”. The last line uses the cumu- 
lant expansion theorem (e^) = ^ and ([/i;A]'^)c((?) denotes the cumulant 

evaluated ai q = qi — q 2 - 

These cumulants can be calculated by expanding the displacement in powers oi 5 l 
and solving eq. (2.22) as in [14]. However, in ID all orders beyond linear order are zero! 
There are two ways of seeing this. The more technical approach is to note that in ID, 
1 + h = det[l + and V = (1 + such that taking the gradient of eq. (2.22), 

noting that = {3/2)H‘^6, yields 


(1 + V,T(g))-i V, fT(g) + 2A4/(g)l = (l - det[l + V,T((z)]-i) . 


(2.28) 


Because in ID det[l + V^T] = 1 + the above equation reduces to 


V„ 


T((?) + 2A4/(g)J = -H^Vg^iq), 


(2.29) 


which is a linear equation and, hence, the linear-order solution is the exact solution of LPT. 
This does not mean that higher order cumulants do not contribute to the true nonlinear 
power spectrum. The expression 1 -|- J = det[l -|- becomes infinite at shell crossing 

and, hence, LPT only applies up to shell crossing (which occurs when (5 ~ 1). 

A more physical way to see that the lowest order solution is the exact solution up until 
shell crossing is from the behavior of the force for an infinite sheet of matter. Since the force 
is directed towards the sheet but independent of the distance from the sheet, prior to shell 
crossing we can replace V(/)(g + 'L) in eq. (2.22) with 'V(j){q), which results in a linear equation 
and so the lowest order solution is exact. 


2.2.1 The Zeldovich approximation 

Linear (or first) order Lagrangian perturbation theory is called “the Zeldovich approxima¬ 
tion” [1, 54, 55]. The Zeldovich approximation is the solution to eq. (2.29) and, hence, is 
exact up to shell crossing. Since SPT also fails at shell crossing (as it makes a fluid approxi¬ 
mation) the exactness of the Zeldovich approximation suggests that in ID LPT should fare 
at least as well as SPT evaluated to any order. In fact, we will show that ID SPT converges 
to the Zeldovich approximation with increasing order. 

The Zeldovich displacement for a particle initially at (Lagrangian) position q is 


'I'za(9) = 


Following eq. (2.25), the overdensity is 


dk jh n / 7 \ 

S' 


^ZA = 


dqe-^^^ - l) , 


-ikq [-ifcTzA(g)]" 

n\ 




n=l 




E 


(2.30) 


=i (2^ 


(2.31) 

(2.32) 

,kn)hiki) ■ ■ - hikn), (2.33) 


^ 2=1 
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where for Fn it follows from eqs. (2.30), (2.32), and (2.33), and from Appendix B.l for 


that [56] 


,kn) = G^r{kir-- ,kn) = -, 




n\ Y\7=iki 


(2.34) 


with k = ki + ■ ■ ■ + kn- We have labelled the Zeldovich kernel as and to identify 
them with the symmetric kernel in SPT (eqs. 2.10 and 2.9). This identification is proven 
in Appendix B.l. Hence, as n —>■ oo the ID SPT solution converges to that in LPT (the 
Zeldovich approximation). The result that SPT converges to LPT is often taken for granted 
in perturbation theory studies (in 3D), as [20] and [57] have shown respectively that to third 
and fourth order in 6l the two theories agree. This result that the expansion of LPT yields 
the same as SPT is nonetheless surprising: The dynamical equations for the two theories are 
quite different (with LPT allowing streams of matter unlike SPT). Nevertheless, the pertur¬ 
bative solutions to these equations are identical. 


All statistical quantities in the ID Zeldovich approximation can be calculated from the 
variance of the differential displacement between two points separated by distance q: 


^^(9) = (['I'za(O) - 'i;'zA(g)]^) = 


dk 2 Piik) 


TT 


k^ 


(1 — cos[A: g]). 


(2.35) 


For example, the Zeldovich approximation power spectrum follows from eq. (2.27) [58-60]: 

PzAik) = j dqe-^'^^{e 

/ OO 

dqe^’^^J2- 


n=l 


,-k^P{q)l2 _ ^ 

(2.36) 

[-A:V(g)]” 

2”re! 

(2.37) 


where we have used that the displacements are Gaussian random and so only the second 
cumulant is nonzero and equal to and the second line simply expands the exponential. 
Since LPT has identical Fn as in SPT, it follows that the re = 0 term in eq. (2.37) yields 
linear theory, the re = 2 term is identical to 1-loop SPT, and so on. With this identification 
between order in and the loop in SPT, once is pre-tabulated eq. (2.37) can be used to 
calculate the SPT contribution to any order in ID with a single integral over q (and indeed, 
for numerical evaluation, a sum of Fast Fourier Transforms). In addition, the correlation 
function is just the Fourier transform of the power spectrum and given by [27] 


iQ-r\\)" 

2cj2(g) 


(2.38) 


Despite its apparent simplicity, because of the q dependence of a the RHS has to be evalu¬ 
ated numerically with rare exception. Equation (2.38) can also be expanded in and its 
derivatives, with each increasing order being one higher loop in SPT (Appendix B.2). 


2.2.2 Common resummation schemes 

In the decades since standard Eulerian and Lagrangian perturbation theory were devised, 
the concentration has predominantly been on schemes to accelerate the convergence of these 
original perturbation theories by resuming terms in the original expansions [17-20]. Since 
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these resummation theories are based on the same dynamical equations as standard pertur¬ 
bation theories, they still rest on the standard expansion converging to a meaningful result - 
one question we aim to address by studying ID dynamics. Their infinite order solutions are 
indeed the same as this limit in SPT and LPT. 

The “improvement” of many of these schemes is in essence to write eq. (2.37) without 
expanding the zero lag term ry = so that in ID 

/ CXD ^ 

(2.39) 

n=l 

where at “linear” order (the n = 1 term) P{k) = exp[—This form appears in the 
Lagrangian theory of Matsubara [20] and an analogous form in Renormalized Perturbation 
Theory [RPT, 18]. In RPT, this zero-lag exponential suppression factor was championed 
as a feature of the theory that it “shields” from small wavelength modes. However, this 
exponential suppression is artificial, occurring because only part of a'^{q) has been kept 
exponentiated; the other part that tends to cancel it has been expanded (note cr^(g) —>■ 0 as 
(y —7- 0). Physically this expansion breaks the Galilean invariance of the theory at every finite 
order: small-scale perturbations do not depend on the properties of large-scale flows at any 
order in SPT or LPT [25, 27, 30, this becomes evident in later expressions]; however, the 
exp[—A:^r/^] term introduces such a dependence. While P{k) = ex.p[—k‘^r]‘^]PL{k) reproduces 
the suppression of the BAG in the correlation function in the concordance cosmology, this 
success is largely a coincidence. For example, for power-law 3D cosmologies with —3 < n^d < 
— 2, the non-exponentiated form leads to the nonsensical result ry = oo and, hence, P{h) = 0 
at finite order even though LPT is still convergent. 

2.3 Comparison of standard perturbation theories 

Here we provide a physical sense for how the aforementioned theories behave before we turn 
to newer effective theories. Since we showed that SPT, LPT, and resummation schemes 
converge to the same result, in a sense these are all the same theory. Figure 2 compares the 
predictions of linear theory, the Zeldovich approximation (equivalent to LPT at any order), 
2 nd spx, and the full nonlinear solution calculated with our Wbody code (section A). 

We show the predicted evolution at a = 0.1, 0.5 and 1 for both a Gaussian perturbation 
with initial standard deviation of 10 Mpc that achieves <5^ = 2 at its peak when a = 1 (top 
panel) as well as our GDM-like cosmology where we have additionally damped the linear 
power spectrum by exp[—5fc^] to eliminate small-scale structures (bottom panel). Especially 
in the bottom panel of figure 2, we see that 2"^^ order SPT performs very poorly relative 
to the full nonlinear solution. Even smoothing the SPT field on ~ 10 Mpc scales - scales 
at which this theory is frequently used - will not cause it to fair nearly as well as LPT. 
Both panels in figure 2 show that the Zeldovich approximation performs excellently in most 
locations. However, in nonlinear locations where shell crossing has occurred the structures 
in the Zeldovich approximation are more extended than in the simulations. The Zeldovich 
approximation also tracks the advection of structures much better than either linear theory 
or 2“*^ order SPT, which is very apparent in a movie of the temporal evolution in the case 
shown in the bottom panel. In such a movie, the linear theory structures are hxed in position, 
whereas the structures in LPT and in the simulation advect substantially across the 100 Mpc 
frame. 

A more quantitative comparison of these theories can be done using the power spectrum. 
Figure 3 shows this calculation for our GDM-like ID model and for various orders in SPT, 
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Figure 2. Top panel: The evolution of a Gaussian perturbation that in linear theory has standard 
deviation of 10 Mpc and an amplitude at peak of = 2 for a = 1. The “nonlinear” curves are 
computed using our particle mesh code, whereas the other curves are the predictions of Eulerian 
linear theory, the Zeldovich approximation, and second order Lagrangian perturbation theory. Bottom 
panel: The same calculations are shown but for a realization of our CDM-like cosmology and where 
we have additionally damped the linear power spectrum by exp[—5fc^] to reduce features. This frame 
shows only 1% of the 10"^ Mpc box used to calculate these fields. 


including the 20-loop power spectrum (i.e. the solution to order As the SPT order 

increases, the solution tracks the Zeldovich approximation prediction to higher and higher 
wavenumbers. Each higher order in the expansion has terms that appear with a higher power 
of k than the previous order, which manifests in the high frequency oscillations that are seen 
at k > 0.1 Mpc“^ in the > 5-loop order calculations. One difference between the ID and 
3D cases is that ID is guaranteed to converge on all scales to the Zeldovich approximation, 
whereas in 3D there is no reason to believe that the series converges to anything finite on 
scales where (5 > 1. However, even though the series converges in ID, it does not converge 
to the result of the A^-body simulation given by the dot-dashed curve in figure 3. Somewhat 
surprisingly, even though the Zeldovich approximation describes the position of structures 
in figure 2 very well, the power spectra differs considerably from the simulation starting at 
k ~ 0.02 Mpc“^. Because the Zeldovich approximation is exact until shell crossing, this bias 
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Figure 3. Top panel: z = 0 matter overdensity power spectrum in our ID CDM-like model calculated 
analytically using linear theory, the Zeldovich approximation (LPT at any order), and SPT to the 
specified order in the overdensity. Note that infinite order SPT yields the Zeldovich approximation 
(which is exact up to shell crossing). The dot-dashed curve is the nonlinear evolution calculated 
from a simulation using this cosmology. None of the calculations converge to the simulation result 
except on scales where linear theory holds. Bottom panel: Shown are the fractional residuals with the 
nonlinear calculation of select perturbation theory curves from the top panel. 


must owe to the deficiencies of the perturbation theory once shell crossing occurs. This shows 
that small errors on ~ 1 Mpc scales propagate to much smaller k, highlighting one of the 
difficulties of a successful perturbation theory. 

2.4 Effective theories of large-scale structure 

We have already noted that the traditional perturbation theory expansions treat nonlinear 
scales as if they were perturbative and treat the cosmological matter field as a perfect fluid. 
The goal of the effective held theory of large-scale structure is to overcome these dehciencies 
by formulating a theory for the dynamics of perturbative (large-scale) modes. These theories 
attempt to write dynamical equations for the smoothed matter overdensity and velocity helds, 
where the smoothing is on large enough scales that perturbation theory can be rigorously 
applied [21, 24, 25, 31, 61, 62]. The effects of the short-range dynamics on the smoothed helds 
is encapsulated in a series of parameters. In the next subsection we detail the effective theory 
approach in an Eulerian context. (See also Appendix C.2 for a description of an approach that 
derives the effective equations based just on symmetries, summarizing [62], although these 
symmetries are also taken into account in the standard approach.) A pseudo-Lagrangian 
effective theory, and a comparison to simulations are treated in ensuing subsections. 
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2.4.1 The (Eulerian) effective field theory of large-scale structure 

The following summary of the effective field theory of large-scale structure (EFTLSS) involves 
several steps, but the final result at 1-loop order is simple and follows from adding a new 
term to the Euler equation that at lowest order in bju is not forbidden by symmetry (and 
physically corresponds to including effective sound and viscosity terms). Let us define the 
smoothing operation on field X as 


Xi = J dxWAix-x')X{x'), (2.40) 

where Wa is a window function with comoving width chosen such that A < A:nl- Here 
/cnl is the “nonlinear” wavenumber below which 5 < 1 such that the smoothed field can be 
treated perturbatively. The EFTLSS as presented in [21, 31, 63] assumes that dynamically 
the universe is a system of collisionless particles. The equations for the dynamics of the 
smoothed fields are 


dJi + ei + x{6m) = -V[u(<5 - 5z)]«, (2.41) 

drOi + nei + i7rGa‘^5i + V{uiVui) = -p’^VVa - Vdr{[ui6 - (iz)]z(l - 6i)} (2.42) 

-V2 {ui[ui5 - 5z)]0 - nV{[ui5 - ,5z)]Kl - 6i)}, 

both valid to 0{6f, A~‘^) for the short-wavelength sensitive terms on the RHS, the order 
required for a 1-loop calculation as demonstrated later. Here 

ta = [ah - piuf] - [[{V<pf]i - (V0z)2] , (2.43) 

pi, ui, au,i, and p are the large-scale density, peculiar velocity, peculiar momentum dispersion, 
and mean matter density. Appendix C.l derives these equations in ID for a set of gravita¬ 
tionally interacting sheets, starting from the Vlasov-Poisson equation. Both eq. (2.41) and 
eq. (2.42) are similar to the continuity and Euler equation used in SPT except (1) the equa¬ 
tions are in terms of the smoothed (or long wavelength) fields and (2) there are extra terms 
that are sensitive to small-scale perturbations. The extra terms in our equations are different 
from earlier attempts, such as [21, 31], in which only the “stress tensor” term, r, appeared 
on the RHS of the Euler equation and no terms on the RHS of the continuity equation. This 
difference only owes to trying to formulate our equations in terms of the Eulerian velocity 
rather than the less physically meaningful velocity that appears in these other papers (see 
Appendix C.l). This reformulation does not have any impact on the final results, but we do 
use this formulation when we estimate the free parameters of this theory from simulations in 
Appendix D.l. 

To solve perturbatively the EFTLSS equations (2.41-2.42) requires calculating the re¬ 
tarded Green’s function of the linear part of these equations (e.g., [18, 31]). First, it is 
simplest to reformulate these equations as a second-order differential equation in just bp. 

— o?T-L^{a)d‘lbi — a {2'H^{a) -|- a'H{a)dH{a)/da) dabi + ^i^Gpa^bi = 

{aUda + n) X{biui) - X{uiXui) 

-P-'vVa + VaUdaOuib - bi)]ibi) - ([u((i - bi)]iui) - Vn[u{b - 50]z(l - Si), (2.44) 

'-V-' 

-P-1V2XA 
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where the final line, which we define as (such an X\ exists), are the terms that do 

not appear in SPT and depend on short-wavelength modes.® The Green’s function, G{a,a'), 
of the linear part of this equation satisfies 

—a?T-L^{a)d‘lG{a, a')—a (2'H‘^{a) + aH{a)d'H{a)/da) daG{a, a)+4:TTGpa‘^G{a, a') = 5^{a—a'). 

(2.45) 

The retarded Green’s function can be written as a linear combination of the growing (oc a) 
and decaying (oc a“®/^) mode solutions with a Heaviside 6^ function enforcing the causality- 
condition G(a, a') = 0 for a < a''J 


G{a, a) 





a 

a' 


(2.46) 


The aim is to be able to solve eq. (2.44) perturbatively. To do so, we expand all new 
small-scale-sensitive terms (which are all not in SPT and hence unique to EFTLSS) in the 
long-wavelength fields, 6i and ui.^ This expansion should depend on all combinations that 
do not break translational invariance and isotropy: 

Xa = Pes + pcldi - -t- ..., (2.47) 

where we have also only kept terms that are lowest order in V"". Here, we use the same nota¬ 
tion as [21] for the coefficients of this expansion - peff, c^, and - because, as shown there, 
there are analogues between these parameters and pressures, sound speeds and viscosities, 
and the Navier-Stokes equation for an imperfect fluid, especially when considering only the 
ta component of Xa- Indeed, we will refer to the effective correction to the power spectrum 
we derive as owing to an effective sound speed. In addition, J is a stochastic component that 
does not correlate with the long-wavelength fields. These parameters, which encapsulate the 
short-wavelength physics that has been “integrated out”, depend on both a and A. 

There is one complication with our derivation so far (see also Appendix G.l). Long- 
wavelength modes depend on the small-wavelength modes at all previous times (e.g., at all 
times along a trajectory; [24] and maybe even on 5™ of these modes). However, the expansion 
of Xa in terms of long wavelength modes (eq. 2.47) assumed locality in time: that the rate 
of change of (5;, ui) at time t only depends on short-wavelength modes at t. Fortunately, at 
1-loop order the final result is equivalent if we treat eq. (2.47) as the fundamental expansion 
as non-locality in time can be absorbed into the time dependence of Cg, c^, and J. These 
complications make the predictions of the previous > 1 loop FFTLSS calculation depend 
on an unknown temporal response kernel [22],® although it is quite possible that future 
formulations will resolve this deficiency [64]. 

®To eliminate ui so the equation is just in terms of Si, we would need to continually plug in w; = J 6i = 
— JldrSi -I- \/{SiUi) -I- V[M((i — 5;)]i] into the RHS of the equation. 

^One can construct a Green’s function from the solutions of the homogeneous equation by taking the 
linear combination of the growing and decaying mode that vanishes at a = a' and then by multiplying this 
combination by 6^^. 

®Up until this expansion, the coarse-grained theory of [62, 63] is identical to the EFTLSS approach of [31]. 
Instead of expanding, this theory measures ta using simulations. 

^Different terms that contribute at 2-loop order have different integrals over this response and so, unlike 
the 1-loop case (where it can be absorbed into a single coefficient), the form of the temporal response affects 
the fc-dependence of the power spectrum. 
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We follow [31] and treat all terms in X\, including as nonlinear terms. The 
motivation for this will be discussed shortly, but boils down to, e.g., serving to cancel the 
UV sensitivity of terms that scale as 5\ and so should be thought of as a third order term. To 
generate the nonlinear solution to the dynamical equations requires integrating G{a, a') times 
the nonlinear terms that appear in eq. (2.44). Since the nonlinear terms depend on unknown 
fields 5i and 9i, a perturbative solution is required. We adopt the standard perturbation 
theory approach, first solving the smoothed linear theory differential equation, yielding 
We then feed this solution into the quadratic terms to solve for the next order solution and 
so on, as described below.^° Note that at lowest order Vui equals —dr log D{a)6L,i and that 
Si equals S^^i, and so up to stochastic and higher order terms X\ is proportional to S^^i. 

When computing the power spectrum with this expansion method, the ETHLSS solu¬ 
tions to order Sf (i.e. 1-loop) include the SPT solutions P 22 and P 13 with 5^ —)• Sl^i. One 
can derive these contributions to .Peftlss^^) taking the nonlinear terms in the eq. (2.44) 
evaluated with the lower order solutions and integrating them over the Green’s function [31]. 
In addition, the other contribution at 1-loop that is unique to EFTLSS is computed in the 
analogous manner and given by [31] 


i2 r ^ f 2 , , drlogD{a') 2 ^ /^ \ j n 

[Sc\i =k I da G{a,a) [Cg{a)-\ -- Cr{a ) \ 5L,i{k,a ). 


It is useful to define 


nia' 




ac = D{a) ^ daG{a,a') |^cf(a') -F ^ D{a'), 


-- i in the Einstein de Sitter if Ctot(aO oc a, 
9 T-L\a)‘^ 


(2.48) 


(2.49) 

(2.50) 


where the approximation comes from requiring that acPii has approximately the same 
dependence as P 13 , as is required to cancel the A dependence. This was the scaling assumed 
in [31], but there will also be a contribution to ctot that has a different time dependence [61]. 

Combining all the terms, the familiar SPT terms and the new terms owing to “sound” 
and “stochasticity”, the 1-loop power spectrum in EFTLSS is 

'^EFTLSs(^) ~ + P22,A + Pl3,A + 2Q;c,Afc^Pll,A + Pj.Aj (2-51) 

where subscript A is shorthand for the replacement P^ —?■ Pl in that terms integrals 

or, in the case of the stochastic term, Pj —?• PjW^. The third and fourth term are the 

--(3) ~ 

only 1-loop terms unique to EFTLSS, and we have used that 2{[dc]i 5l,i) = 2ac,APii,A- In 
contrast to previous EFTLSS studies, we have defined ac to be negative for positive ctot, in 
part because in most ID cosmologies we consider ctot is negative. 

EFTLSS assumes that modes with /c' S> /c do not influence -Peftlss^^) ™ 
except through the value of 2ac- It is easily verified that the new terms can effectively absorb 

^°There is an inconsistency in this way solving the solution in that at each order in 5l,i the solution should 
be smoothed once by our filter, but the filter enters in a more complicated manner [24]. This complication is 
ultimately resolved with our limiting procedure. 

^^The authors of [21] argued that the A = oo sound speed, Cs, is always positive and generally will dominate 
over the contribution from c„ in 3D. However, their argument for the positivity of Cs relied on the nonlinear 
power spectrum growing faster than the linear one, which we find to be never the case in ID. 
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the contributions that SPT receives from k' ^ k. The EFTLSS “sound” term is proportional 
to k‘^PL{k), which is precisely the /c-dependence needed to cancel the UV contribution in P 13 
(i.e. to cancel the A dependence). Also, modes with k' ^ k contribute a term in P 22 that 
scales as k^ as k ^ 0, just like the stochastic term Pj (see appendix B.3). However, prior 
formulations of EFTLSS drop any term that has a weaker scaling in /c as fc —>■ 0 than higher 
loop terms for consistency. For the CDM case, this results in dropping the stochastic term, 
Pj, even at 2-loop order [22]. This may be surprising because halo models for large-scale 
structure explain nonlinear evolution with just a stochastic term (scaling with the abundance 
of halos) that goes as k^ [33-35]. However, in power-law cosmologies one can show that the 
higher order contributions to the stochastic term can be expanded as (A:/A;nl)^'*~^* with order 
unity coefficients for integer i > 0 [61]. Thus, these terms are subdominant at /c <C /cnl- In 
Appendix D, we show that the stochastic term is also subdominant to our 1-loop corrections 
in our CDM-like cosmology on scales where 1-loop perturbation theory is successful. 

Once the stochastic term is dropped and ac is determined, so that all terms in eq. (2.51) 
are calculable, the story of EFTLSS is not over. It would be problematic for the solution to 
depend on A on any scale at which the calculation is used. However, for A that correspond 
to plausible the predictions at mildly nonlinear scales of our low order EFTLSS (as 
formulated above) do tend to depend on A. (For A“^ = lOMpc, [rf]^ is suppressed by 20% 
relative to [r/^]oo while P 22 is less sensitive to A.) To remove this A-dependence requires 
including terms to higher order in fe/A, working at lower k, or increasing A. The approach 
taken in the EFTLSS literature to avoid this A-dependence is to take the limit A —)> 00 in the 
standard integrals and for the EFTLSS parameters (in our case ctot)- While this is the limit 
of a non-perturbative field that the EFTLSS equations were derived to avoid, the assertion 
of EFTLSS is that this limit makes more sense to apply at the final step. (The mechanics 
of how this limit is applied are apparent in the power-law cases discussed in section 3 and 
suggest that this limiting procedure is not unreasonable.) The A —)• 00 limit can be done by 
extrapolating measurements smoothed on smaller and smaller A or by using the coefficients 
in the A = 00 expression coefficients that provide the best fit to simulations. For our ID case, 
this extrapolation is quite simple because the sound term must add to Hi 3 ,a = —^^[??^]aT’l,a 
in a manner that sums to something A-independent, requiring 

2ac,oo = 2ac,A + ([^^]a ~ ['^^]oo) • (2.52) 

In what follows, we drop the 00 subscript and write 2 q;c,oo as 2 q;c. 

In the following subsections, we develop a pseudo-Lagrangian effective theory (Sec¬ 
tion 2.4.2) and test the predictions of EFTLSS on the nonlinear power spectrum from simu¬ 
lations (Section 2.4.3). (Additionally, in Appendix D.l, we use our ID simulations to estimate 
2ac and hnd consistency with the power-spectrum matching technique used in the main text.) 
Foreshadowing, we find that the EFTLSS is remarkably successful. 

2.4.2 A Lagrangian effective formulation 

As with SPT, a deficiency of Eulerian EFTLSS is that this theory expands the (smoothed) 
matter displacement, af, in density [25]. Since the infrared displacement is close to being 
non-perturbative at the BAO scale despite coming from modes where linear theory is highly 
applicable (which we demonstrate in detail later in section 4), this leads to a large error in the 
predictions at this important scale. This motivates again constructing a “Lagrangian-space” 
effective theory, which we will refer to as LEFTLSS, that resums the contribution to matter 
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displacements from large scales, as done in [25, which related it to a theory of extended 
objects]. We describe a simpler but related Lagrangian approach here. 

Let us image a long-wavelength and short-wavelength decomposition of the Lagrangian 
displacement such that 'h = 'I'; -|- 'Ls- (Somewhat confusingly, our short wavelength term 'L* 
includes equally-long wavelengths as 'h; but sourced by the coupling of small-scale, k > A, 
modes.) We can write the short-wavelength term as an expansion in the long-wavelength 
fields plus a stochastic (long-wavelength uncorrelated) term 

~ 2ac,A^6i + V J, (2.53) 


where the above expansion keeps the lowest order terms that are local in the long-wavelength 
fields (and, importantly, the order required to generate the 1-loop EFTLSS counter terms). 
The coefficient of the former term, which encodes the long-wavelength response of couplings 
with small-scale modes, has been named in a familiar manner because of a later identification. 
In addition, we have dropped terms that are higher order in 6i or its derivatives, and we 
omit a term as at the lowest order of interest it is degenerate with Vh/. We can now 
take eq. (2.36) for the power spectrum in terms of the displacement, which depends on the 
cumulants of 'l'(gi) — T(( 72 )- At our working order in 5i, it only makes sense to consider the 
second cumulant (remembering that in ID LPT only the second cumulant contributes): 


([T(gi)-T(g2)]') = ([^K0)-^z(9)]')+2[2ae,A+4«2^V2][ei,A(0)-eL,A(g)]+2V2[O(0)-O(g)], 

(2.54) 

where again q = qi — q 2 and is the correlation function of J. Without worrying about 
whether UV sensitive terms are countered by other terms (which we consider shortly), let 
us take the limit A —)• oo. As with this step in EFTLSS, our small scale terms do not go to 
zero, such that the A —)• oo power spectrum is 


P{k) = I dqe {e ^^'^left/2 _ ^ 


where CTleft ^ - ^iQ 2 )Y 


A^OO 


(2.55) 

(2.56) 


The above includes all orders in the displacement and the 1-loop order effective theory terms. 
Since the aim of LEFTLSS is to keep the linear IR RMS displacement to infinite order, but 
then to still have consistent counter terms, it makes sense to expand the exponential in all 
UV sensitive terms but to keep the IR RMS displacement. To proceed, we define the variance 
in the displacement contributed by modes that are less than or greater than fc* as cr^ and , 
where /c* is a free parameter to which the results are hopefully relatively insensitive. We now 
want to expand all the UV sensitive terms so that they only enter at their lowest (1-loop) 
order (i.e., keeping Uleft subcomponent of o'left)’ order at which they are 

consistently countered by our terms: 


pl—loop 
-^LEFT 


(k) = Pl^{k) + k^ 




/ {P22+Pl3}> 

—it^/2-|- k'^a^/S -|-2a;c^L + 


V 



where P'^^{k) is the LPT power spectrum calculated using modes only with k < /c*, we 
have dropped the term because it is suppressed by derivatives, and M is a constant that 
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Figure 4. Comparison of the power spectrum in 1-loop EFTLSS to the power-spectrum in the sim¬ 
ulation, in Eulerian linear theory, and in 1-loop SPT. The left panel shows z = 1, and the right panel 
shows z = 0, both in our CDM-like ID cosmology. The red horizontal lines span the wavenumbers 
over which 2ac is estimated, although it makes little difference if a factor of two smaller wavenumbers 
are used (as apparent from the flatness of the residuals). In the right panel, we also show the Zeldovich 
approximation power spectrum and the predictions of 1-loop LEFTLSS (discussed in section 2.4.2) 
where 2ac is fit separately again as a free parameter. 


encapsulates ^l(O) and other terms that renormalize it^^; it is suppressed at nonzero lag by a 
factor of P'^^{k) compared to other terms. Notice that all UV divergences are now properly 
cancelled. This formula is similar to the expressions derived in [25]; see also [32, 65]. The 
limit /c* —)■ 0 with ^ = 0 of yields T’gp:^Lss ■ 

In addition, we can calculate the correlation function by Fourier transforming eq. (2.57) 
in the case where A:* is not a function of k: 


^1—loop 
SLEFT 


(r) 


tLPT 

S< 


+ 


(r) + 


\/2 


dq e 


(r-q) 


0-2 - (r - qf 


vrcr; 




3it< — 6(r — q)^a^ + — q) 




(~< 7>/2 + ‘^Oic^L + 

(2.58) 


We subsequently set A to zero, but it does not have to be, and on scales where the stochastic 
term is white then = 0. We note that in the fc* = oo case - which we later find produces 
nearly identical predictions to more motivated A:* choices - , (r) and Fleft^(^) 

particularly simple forms, the sum of the Zeldovich approximation prediction plus a Gaussian 
integral over 2a^L + V^^j. 


2.4.3 Comparison of effective theories with standard approach 

Figure 4 shows P{k) from the simulation, linear theory, 1-loop SPT, and 1-loop Eulerian 
effective theory (EFTLSS) for our CDM-like case at z = 0 (right panel) and z = 1 (left 

Higher order effective terms in the displacement expansion includes terms that when expanded scale in 
the same way. 
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panel). In the right panel only, we also show the Zeldovich approximation power spectrum 
and the predictions of 1-loop LEFTLSS. To determine 2ac, which is needed to compute the 
EFTLSS curves, we have fit for it using the EFTLSS formula over the range shown by the 
red horizontal line, rather than measure it from the small-scale behavior of the simulations. 
The best-fit 2 q;c we find to be 39 Mpc^ for z = 0 and to be 6 Mpc^ at z = 1, which we 
note is somewhat different than the scaling often assumed. For both redshifts, EFTLSS 
dramatically improves the range of scales at which 1-loop perturbation theory is 1% accurate: 
SPT achieves this accuracy over k < 0.05 Mpc“^, whereas EFTLSS achieve it over k < 
0.1 Mpc“^ for the z = 1 case. EFTLSS shows a comparable relative improvement over SPT 
at z = 0. Also, note the EFTLSS residuals are relatively flat over the scales where there is 
improved agreement: This means that the k‘^P{k) scaling is correctly capturing the slope of 
the excess term that 1-loop SPT is missing. It may be surprising that Zeldovich fares much 
worse than EFTLSS, as it is exact up to shell crossing (see right panel in figure 4). This 
improvement implies that EFTLSS, where the primary difference is that it fits for a term 
that scales k‘^P{k), is doing a better job at capturing the impact of regions where shells have 
crossed on the large-scale power. 

The 1-loop Lagrangian effective theory (LEFTLSS) curves in the right panel of figure 4 
are computed by dropping the stochastic contribution (as also done in EFTLSS) and by 
either taking A:* = 0.1 Mpc“^ (using a Gaussian cutoff of the form exp[—A;^//c^]) or fc* = oo. 
The former fe* corresponds to roughly the wavenumber where the Zeldovich approximation 
power spectrum errs at 10%, and /c* = oo corresponds to the advection owing to linear modes 
over all wavenumbers being resumed. However, we find that both values of A;* yield nearly 
the same prediction for our CDM-like case. For the LEFTLSS curves we again fit for 2ac, 
with the fit preferring a different value than in EFTLSS of 2 q;c = 25 Mpc^. That LEFTLSS 
is best fit with a different value is not surprising because the small-scale terms that 2ac is 
countering are somewhat different between the Eulerian and Lagrangian formulations. 

If you squint, LEFTLSS appears to take out the BAG wiggles that are tentatively 
present in the EFTLSS residuals in figure 4, as others have suggested a Lagrangian theory 
should [32, see section 4]. In addition, LEFTLSS fares as well at predicting the power spec¬ 
trum as EFTLSS. We consider the correlation function in section 4, where the improvements 
of LEFTLSS over Eulerian EFTLSS are even more notable. 

3 Power-law models 

Thus far we have considered our CDM-like model. We now consider power-law power- 
spectrum models where the ID linear-theory power spectrum is 

kPLik)/{27T) = (A:/A:nl)”+^ (3.1) 

for k > 0 and the symmetric form for k < 0. During the Einstein-de Sitter phase, the 
evolution of such scale-free initial conditions is self-similar and so it is sufficient to consider 
a single epoch. The panels in figure 5 shows the predictions for power-law models (all with 
^NL = 1 Mpc“^) with index n = 0, n = 0.5, n = l and n = 2, from left to right respectively. 
The curves correspond to linear theory, the Zeldovich approximation (LPT at all orders), 
and the fully nonlinear calculation, all computed by taking the power spectrum of modes in 
the simulation volume (i.e., using the same initial random numbers as the simulation curve; 
^Nyqmst = 30 Mpc“^). We have tested that the nonlinear power spectrum from the simulation 
(calculated using the fiducial specifications that were chosen with the CDM-like cosmology in 
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n=0 n = 0.5 n = l n = 2 



Figure 5. Comparison of perturbation theories and simulations for models with a power-law linear- 
theory power spectrum of the form fcPi,(fc)/(27r) = (fc//cNL)"^^ for k > 0, with /cnl = 1 Mpc“^. The 
dashed line shows Pl, and the dashed-dotted line shows the result of the fV-body simulation. The 
Zeldovich approximation power is shown as the thick black solid line (computed from the modes in 
the simulation with fcNyquist = 30 Mpc“^), while that of 1-loop SPT is shown as the dotted line (using 
a top hat Wa with A = 10). The predictions of EFTLSS are shown as the green solid line, where 
EFTLSS is the same as SPT except that we have fit an extra term that scales as k‘^PL{k) ~ k'^ in 
the left panel, k^PL{k) ~ k'^'^ in the left-middle panel, k^PL{k) ^ k^ in the right-middle panel, and 
Pj k^Phik) ^ k^ in the right panel. The grey band in the rightmost case is to indicate that this 
particular numerical simulation is not fully converged with the estimated error being roughly half the 
vertical width of this band. 


mind) is converged for the three leftmost cases (appendix A), but, even using much smaller 
time steps and higher initialization redshifts, we were only able to achieve convergence at 
the 2% level for the rightmost n = 2 case.^^ Figure 5 also shows the predictions for P{k) of 
an analytic calculation using 1-loop SPT with a cutoff at A = 10 /snl (a cutoff is required 
to yield a finite answer in the rightmost two panels) and of EFTLSS. The 2 q;c parameter 
in EFTLSS has been adjusted, by hand, to provide a good fit to the A^-body results, and 
we have ignored the stochastic term as is justified for these cosmologies.^^ We do not show 
curves for LEFTLSS as this model makes less sense to apply to these power-law cases where 
the large-scale displacement is not necessarily the effect that is the least perturbative. 

Recall that both SPT and LPT can be expressed as expansions in cr^ (q). For power-law^^ 


put some effort into simulating for n = 3, but were not able to achieve even a modicum of convergence. 
Indeed, it is somewhat surprising that the n = 1 and n = 2 converge to a steady state (which we take to be 
the self-similar solution) even though the Zeldovich approximation power spectrum that is used to initialize 
the simulation is formally infinite for n > 1 (although not our simulated case that is cutoff aX k — IO/cnl). 
Because n < — 1 cosmologies have large-scales collapse before small scales such that perturbation theory is 
not valid, the range of power-law cases we select essentially maximizes what is possible to simulate. 

^"^For n = 2 the stochastic term is potentially important, but it has exactly the same form as the effective 
sound speed term. 

^®One can find expressions for P{k) for power-law models for the full 3D case in the Zeldovich approximation 
in [66] and within SPT in [61]. 
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models cr‘^{q) converges for — 1 < n < 1, yielding for our normalization convention 


cr^(g) = ^ 9 ^ ’^r[n-l]sin 


(n + 2)7r 


(3.2) 


Infinite order SPT and LPT are thus convergent over the same range in n. That ci^ diverges 
for the rightmost two panels in figure 5 perhaps explains why the Zeldovich approximation 
fares poorly when calculated from the modes within the box. However, the Zeldovich approx¬ 
imation also deviates significantly from P{k) in the simulation for the n = 0 and n = 1/2 
cases, where cr^(g) is convergent. Even Eulerian linear theory matches most of the simu¬ 
lations to higher k than the Zeldovich approximation. This again demonstrates that shell 
crossing plays an important role in the nonlinear evolution of the power spectrum, noting 
the exactness of the Zeldovich approximation prior to shell crossing (or of SPT to infinite 
order). 

One convenient aspect of power-law models is that there are simple analytic forms for 
the 1-loop power spectrum in all the power-law cases considered in figure 5. The 1-loop 
EETLSS predictions assuming a top hat W{^{k) with a sharp cutoff at A are 


pl—loop 
-^EFTLSS 


{k) 


27r -f dvr A: (^log ^ - 27rP ^ + ^irk'^ac^K + Pj,h if n = 0, 

27rA:V2 _ ^ 2 Trk ^/2 - A - if n = 1/2, 

2-Kk + 27Tk^ (log [(1)^(1 - i)] + 2ac,A) + Pj,A if n = 1, 

2Trk‘^ + 2'Kk^ [2ac,A — A] — nk^ + Pj,A if n = 2, 

(3.3) 


setting /cnl = 1 (with SPT being the limit A —)• oo and Pj^\ = 0). Only in the n = 0 and 
n = 1/2 cases are the SPT limits of these equations convergent, reducing to = 27r 

(no correction over Pl) and = 27rfc^/^ — SttA;^, respectively. In both of the n = 0 and 

n = 1/2 cases, the effective terms act to cancel the lowest order in k/h dependences, and 
also in both of these cases EETLSS predicts power-law scalings in k than are not present in 
A —)■ oo SPT. [SPT gives terms that scale as while EETLSS adds a term going as kP~^‘^ 
with an unknown coefficient and Pj, which has a leading-order behavior of as A: —)• 0.] In 
the n = 1 and n = 2 cases, the new EETLSS terms act to exactly cancel bone fide divergences 
as A —)• oo.^® In the n = 2 case, both the 2ack^PL and Pj^\ have the same A: —>• 0 scaling (A:'^) 
and so both act to cancel the SPT divergence that scales as A;^A. At higher n, which we find 
are not feasible to simulate, Pj^a rather than 2ack‘^PL cancels the lowest order divergence in 


SPT. 


In all power-law cases considered in figure 5, EETLSS agrees much better with the 
simulations than with LPT (using the modes in the box) or with SPT (cutoff with A = 10), 
providing a factor of ~ 3 enhancement in the scales beyond where linear theory is percent- 
level accurate (or the other theories for that matter). This improvement does not appear 
to be a coincidence because of the additional free parameter. In the two cases where SPT 
yields a finite prediction, EETLSS is still a significant improvement over SPT. While an 
improvement over 1-loop SPT is guaranteed because EETLSS has a free parameter, we see 
that these additional terms are at least as important as the perturbation theory generated 
ones for n > 0, allowing us to fit the A^-body power spectra at the percent level to about 


^®The sound term also counters the A dependence for the logarithmic term noting that this term is 
21og(fc/A). 
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Figure 6. The top panel shows the 1-loop SPT, Zeldovich approximation/LPT, and 1-loop LEFTLSS 
predictions for the correlation function (^), and the first three labeled curves in the bottom panel 
correspond to the 1-loop terms that contribute additively to ^ (using eq. 4.1), with the smaller terms 
multiplied by a factor of ten. The other curve in the bottom panel is ten times the difference between 
LEFTLSS and LPT. Note that the full nonlinear solution agrees remarkably well with the Zeldovich 
approximation (figure 1). These calculations are for our ID CDM-like cosmology and ^ = 1, an epoch 
that shows comparable nonlinear evolution to these functions at z = 0 in 3D. The bottom panel shows 
that the RMS displacement term, is larger than the other 1-loop corrections by more than a 

factor of ten. 


0.3 A^nl or so (and, actually, this is about the fraction of A:nl we fonnd for our CDM-like case). 
Fitting instead for power-law terms of the form fc™ with integer (or half integer) powers in 
m and that are different than the EFTLSS terms does not work nearly as well, a result that 
can be gleaned by noting that the slope of the EFTLSS term is exactly what is required 
to correct the SPT prediction. Pajer &: Zaldarriaga [61] presented a similar comparison of 
EFTLSS against power-law models, testing it for 3D cosmologies with n = —1.5 and — 1 
using data from published simulations. They fonnd < 5% residuals, and, likely because of 
the accuracy of the simulations they employed, their residuals are not as striking in their 
flatness as those shown here. 


4 Understanding the 1-loop behavior of these theories 

Here we attempt to understand the 1-loop predictions of these different theories and how 
well this order captnres the actual nonlinear evolution, again specializing to the ID case and 
returning to CDM-like cosmologies. We show that this simplified case has bearing on many 
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of the topics discussed recently in the perturbation theory literature (such as shifts in the 
BAO peak, the comparison of LPT to SPT, and the [dis]advantages of a Fourier space over 
a configuration space analysis). 

Let us first consider the predictions of these theories in configuration space. Configu¬ 
ration space is the most convenient space to compare with observations of galaxy clustering. 
The 1-loop correlation function is simply the Fourier transform of the 1-loop expansions for 
P{k) (eq. 2.13; appendix B.2): 


>1—loop 
•^SPT 


(x) 


growth 


dilation 


^l{x) + ^^1{x) + 4^'^{x) dx^L{x) + 


RMS displ. 



(4.1) 


where primes denote derivatives with respect to the argument. For 1-loop EFTLSS, we would 
replace a‘^{x) with cr‘^g{x) = (t(x)^ — 2q:c in eq. (4.1).^^ The 2ac term is a 15% correction to 
in our CDM-like cosmology, and we shall ignore this correction in most of what follows 
as this difference has little impact on our discussion. 

Each of the terms in eq. (4.1) can be thought of as the lowest order expansion with 
respect to a certain effect. The first term rightward of the linear term owes to the volumetric 
average of how regions grow relative to linear theory, the second term encodes the dilation 
that large-scale fluctuations impart on smaller separation correlations,^® and the final term 
is the effect of the RMS matter displacement that owes primarily to modes with wavelength 
> X. The bottom panel in figure 6 shows the contribution of each of these terms to a 
quantity that emphasizes the BAO scale. The solid curve in the bottom panel represents the 
dominant correction owing to the RMS displacement, (T^^^/2. The other terms that appear in 
eq. (4.1) [all multiplied by a factor of ten in figure 6] result in much smaller corrections at the 
BAO scale and indeed even smaller corrections to the underlying continuum at x > lOMpc 
for the CDM-like case at hand (all 1-loop corrections scale with time in the same way). This 
result is not surprising when looking at these terms in detail. In CDM, the RMS displacement 
at the BAO scale is iT(lOOMpc) ~ lOMpc, comparable to the width of the BAO feature and 
thus to I/y^^. Hence, the RMS displacement term, (T^^^/2, is an 0(1) correction. The 
growth term, is small because |^L(ic)| rsj 10 ® at X r-u lO^Mpc and z = 1 (see figure 7). 
The dilation term, Ax^i{6‘l)x, is also small because {6‘l)x = x~^ dx^i 10 ® and 
10 at the BAO peak. The same result that the RMS matter displacement is the 

main source of evolution of the BAO feature in the correlation function applies to the full 
3D CDM case [30, 32]. 

can write a similar expression for S{x) to second and third order: 


^(2) 


^(3) 


growth 
^ + 


growth 


(V<5l) [ 

J X 


dx 8 L 

2 ^^^ order shift 


growth —shift 


Si +2 ^"Sl 


/ POO \ ^ POO 

ij dx 5 l\ +ZSlV5l J dx Sl ■ 


(4.2) 


(4.3) 


Somewhat confusingly, both shift and growth terms in 5 can mix to generate what we call a, e.g., growth term 
in ^{x) (eq. 4.1). 

^*In overdense regions, the large-scale overdensity acts like a closed universe remapping r to slightly smaller 
scale. This effect does not cancel when averaging over underdense and overdense regions because there is 
more growth in overdense regions, such that they are weighted more heavily than underdense ones (as we 
show shortly). 
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Despite its smallness, the dilation term has received some attention because it can shift 
the peak position of the baryon acoustic scale [26, 40, 67], whereas the other nonlinear terms 
do not result in a significant shift because they are quadratic in derivatives of or, in the 
case of CT, have a broadband response. To compute the BAO peak shift from the dilation term 
we note that .^^(x + ex) ss + ex^'i{x). Thus, the dilation term results in a displacement 
in the BAO peak in the correlation function of e ~ 4(5|^(x ~ lOOMpc)). The ID shift using 
1-loop SPT is 3.2 times larger than the shift of 68/63 (j£(r ~ lOOMpc)) found in [40] for 
the full 3D case. The authors of [40] showed that the 3D case can be understood from the 
following factors; that the local correlation function is enhanced at lowest order by a large- 
scale overdensity, 5 lS) by (1 -|- ss 1 -|- section 2.1, and the correlation function 

is shifted by large-scale overdensities as « ^l([ 1 + <^Ls/3]r) owing to comoving scales being 
contracted by the large-scale overdensity. Thus, in 3D, expanding around r, 

?(r) ~ (^CUr) + • V^L(r)^ ^ Ri ^^(r) -h ^(<^ls) ^ 

Repeating the same exercise in ID, the correlation function is modulated locally by (1-|-2(5ls)^ 
(section 2.1) and the ID shift is |l([1 + <5ls] 3;). Combining these effects in the same manner 
yields 4 ((5^g) exactly our ID dilation term. 

Even at z = 1 in our fiducial CDM-like cosmology, 1-loop SPT fares poorly at describ¬ 
ing the shape of the BAO in the correlation function (top panel; figure 6).^® The lack of 
success owes to 1-loop order dropping terms that are higher order expansions in the RMS 
displacement, of the form iT^’"[V^]™'^i(x). We have checked that the term cr'^^^"(r)/8, which 
enters at 2-loop order, is not substantially smaller than the dominant nonlinear term in the 
correlation function a^^'^'{r)(2 at z = 1, and it grows in importance relative to lower loop 
terms with time. This demonstrates that SPT is only going to be weakly convergent in 
the correlation function at BAO scales. In contrast, LPT includes all orders of these RMS- 
displacement terms (in both ID and 3D) and hence should fair much better. This inclusion 
was also the motivation for LEFTLSS. Figure 6 shows the 1-loop LEFTLSS calculation us¬ 
ing our best-fit EFTLSS value of a = 6 Mpc^ and = oo (section 2.4.3; although we find 
using /c* = 0.1 Mpc“^ is essentially indistinguishable).^® The LEFTLSS prediction is nearly 
identical to the very successful Zeldovich approximation one, further confirming that the 
UV-sensitive corrections are small at BAO scales. 

Now we turn to Fourier space and the 1-loop corrections to the matter power spectrum. 
The density, dilation, and RMS displacement terms that constitute the 1-loop corrections in 
the correlation function (eq. 4.1) each have a power spectrum that is generally much greater 
in absolute magnitude than the linear-theory power spectrum at all wavenumbers; there are 
large cancelations when they are summed into the 1-loop SPT power spectrum. However, we 
can rewrite the terms in a manner that removes the terms that cancel, decomposing 
as a term that mixes the different effects and a term that derives purely from the RMS 


^®Even though 1-loop is not terribly successful at capturing the nonlinear evolution at BAO scales in our 
ID CDM-like example, we note that the failure is even worse in 3D CDM where the 1-loop prediction in the 
correlation function is infinite! 

^°Really we should fit for this parameter in LEFTLSS, but this matters little here. At 2 ; = 0, such fitting 
resulted in a 30% smaller a^. 
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Figure 7. Decomposition of our CDM-like ID power spectrum at 2 ; = 1. The labels distinguish 
terms that owe to the RMS displacement from “mixed” terms that mostly arise from other effects 
(see eq. 4.4). Also shown is the absolute difference between the 2-loop and 1-loop power spectrum and 
the same but between the 20-loop and the 1-loop, both multiplied by 5 (labelled 5 x |2-loop| and 5x |20- 
loop|). The RMS displacement effect is comparable to other 1-loop terms for the first and second 
BAO peaks in the power spectrum, in contrast to the correlation function, and the SPT expansion 
also converges more quickly at these scales than at the BAO peak in the correlation function. 


displacement: 


/ oo jy 

^1 


mixed 


RMS displ. 


2^PL{k')PL{k - k') + ^PL{k')PL{k - k') - kWP{k)5^{k') 


(4.4) 

where we could replace rf with rf — 2ac,oo to yield EFTLSS, but we shall neglect this small 
change in what follows. Figure 7 plots the mixed and RMS displacement terms at 2 ; = 1 in 
our CDM-like case. The mixed term is larger than our problematic term for the first two 
BAO peaks (the peaks where SPT is used in 3D), and unlike the RMS term does not show 
acoustic oscillations. Also shown are the contribution of higher loops. In particular, the 
yellow dashed curve is the absolute difference between the 2-loop and 1-loop power spectrum 
multiplied by 5, and note that the terms at each loop order grow as and, hence, 

would be more important at z = 0. The black dashed curve is the same but between the 
20-loop power spectrum and the 1-loop. The beyond 1-loop contributions are smaller at the 
first couple peaks (and including 19 more loops makes little difference at A: < 0.13 Mpc“^ 
for z = 1), suggesting that Eulerian formulation of perturbation theory (SPT and EFTLSS) 
yields a convergent solution for the first ~ 2 peaks in the case at hand. (Of course for ID 
SPT, the convergence is to the Zeldovich approximation solution and not the true nonlinear 
solution.) The lack of convergence in the correlation function at the BAO peak - discussed in 
the first half of this section - owes to contributions from shorter wavelengths at which SPT 
is not trusted in Fourier space. This ID picture differs somewhat from that in [30, albeit 
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in 3D], who argued that the agreement at BAO scales between SPT and the full nonlinear 
solution is spurious, appearing better than it actually is because the continuum power (for 
which SPT is more justified) is comparable to the power in BAO. Lastly, note that the 
contribution from > 1-loop terms show significant oscillations around the BAO, arising from 
the RMS displacement expansion. As in the correlation function, in the power spectrum the 
displacement terms are the ones that converge slowest. 

5 Conclusions 

This paper discussed cosmological perturbation theory in the specialized case of ID dynamics 
(e.g. gravity between expanding sheets). In ID, linear-order Lagrangian perturbation theory 
(LPT) is exact up to shell crossing and all nonlinear terms in LPT are zero. We further 
showed that infinite-order standard Eulerian perturbation theory converges to linear-order 
Lagrangian perturbation theory (and is easily calculable for any given order). Unfortunately, 
for a wide variety of initial conditions, we found that these standard perturbation theories 
do not converge to the nonlinear solution at any mildly nonlinear scale, indicating that 
dynamical equations that allow for shell crossing are needed. 

In this spirit, we considered the recently developed effective field theory of large-scale 
structure (EFTLSS), both the standard Eulerian space formulation and a Lagrangian space 
derivative. This theory attempts to more rigorously formulate perturbation theory by using 
dynamical equations that depend only on smoothed fields whose RMS is < 1. When the 
effective sound speed that appears in EFTLSS was treated as a free parameter, we found that 
at I-loop order EFTLSS extended the range over which perturbation theory is percent-level 
accurate by a factor of two in wavenumber in a CDM-like cosmology in which the overdensity 
has the same variance per k as the concordance cosmology (consistent with the findings of 3D 
studies; e.g. [22, 31]). We also tested EFTLSS in a large range of power-law ID cosmologies, 
where we found its improvements over previous theories were even more dramatic. The 
new terms that EFTLSS predicts based on symmetry are exactly those required to correct 
previous theories. 

Specializing to ID also allowed us to break up the lowest order nonlinear contributions 
to the matter power spectrum into components that correspond to three distinct physical 
effects: the RMS displacement of particles, the effect of an overdensity on its own growth, and 
the modulation of smaller structures by large-scale ones. At the BAO scale in the correlation 
function for a CDM-like case, this allowed us to show that the primary source of nonlinear 
evolution is the large-scale RMS displacement that LPT captures exactly but that are treated 
perturbatively in SPT and Eulerian EFTLSS. In a CDM-like cosmology, we showed that the 
1-loop expressions for these theories provides a poor description for the correlation function 
at the BAO peak owing to the largeness of the RMS displacement, which the 1-loop theory 
has expanded in despite the expansion converging slowly at such scales. This echoes similar 
findings in 3D, but we highlight that the analytics in ID are substantially simplified. In 
contrast to the case of BAO-scales in the correlation function, in the power spectrum the 
RMS displacement term is comparable to other nonlinear terms for the first few acoustic 
peaks. We also discussed beat coupling and the related shift of the BAO peak by large-scale 
couplings in our ID case. 

Of course, the ultimate aim of perturbation theory studies is a rigorous understand¬ 
ing of these theories in 3D. However, our investigation of the dynamics specializing to ID 
has several advantages, including (1) much smaller simulations are required to estimate the 
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nonlinear solution to the required precision, (2) the solutions for standard perturbation the¬ 
ories can be computed to arbitrary order, and (3) the derivations are considerable simplified. 
Furthermore, the analogies between effects in this highly symmetric case and the full 3D 
one are truly striking, making the ID case an excellent means of understanding the complex 
phenomena involved in the formation of large-scale structure. 
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A Particle-mesh code and convergence 

While ID gravity can in principle be integrated exactly as the forces are determined simply 
by the ordering of mass elements [68, 69], we have opted for the simple and fast particle 
mesh (PM) method. The PM method has been used extensively for cosmological A-body 
simulations in part because force computations boil down to Fast Fourier Transforms [70, 71]. 
In the PM method, forces are computed on a grid, and these forces are used to move particles. 
We have written a ID PM code, as also done in [72, 73], that employs the 2°*^ order-accurate 
leap-frog integrator with time steps that are uniform in ^/a, with this choice chosen by trial 
and error from amongst steps with the parameterization for real p. Our code uses CIC 
interpolation to both project the mass elements onto a grid and to evaluate forces from a 
grid of V(/>. In addition, we also smooth the force by an additional CIC kernel such that 
the gravitational potential is damped in Fourier space by sin^(/cAx/2)/(/cAx/2)^ to reduce 
particle noise on the cell scale, where Ax is the size of a mesh element, as motivated in [70]. 
The initial conditions use the Zeldovich approximation displacements. Because the Zeldovich 
approximation is exact until shell crossing in ID, this allows us to start our ID simulations 
at later times than typically done in 3D simulations. Our PM calculations start at z = 10 
unless otherwise specified. 

There are a couple confusing conceptual issues for ID cosmological simulations. First, 
what does it mean to compute forces from an infinite series of sheets (or, in our case, a 
periodic domain with N sheets)? The force on any sheet is proportional to the number of 
sheets to the left minus the number of sheets to the right (both infinite). In a periodic box 
of length L, the force at x is the number of sheets from (x, L) minus that from (x, 0). Thus, 
every sheet in the box feels a force except the center sheet. At first glance this violates 
translational invariance, but note that the Hubble friction exactly counters this gravitational 
force. If sheets have not crossed, how do density structures develop as there is no net force? 
The answer is that sheets start off with different peculiar velocities (as determined by the 6l) 
and, hence, there will be convergences and divergences that affect 6. A final point of clarity 
is that the sheets must have constant surface density in comoving space in order to conserve 
mass. 

To simulate the dynamics of 10® sheets with 10^ PM elements from ^ = 10 to z = 0 
with our code requires about an hour on one CPU. We have parallelized our code over shared 
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- Fiducial; PM=10®“, N=10®®, Box=10’®Mpc 

- High Res; PM=10® ®, N=10®Box=10'®Mpc 

- Large Box; PM=10®’, N=10®’, Box=10”Mpc 

- High N; PM=10®®, N=10®‘', Box=10’"Mpc 



k [Mpc ^ ] 



Figure 8. Convergence tests varying the PM grid size, the number of sheets, the box size, the 
initialization redshift, and the time step for the z = 0 CDM-like power spectrum. Models with fixed 
box size have the same random numbers for all overlapping modes. For all cases we achieve < 1% 
precision. 


memory, which has enabled us to run simulations that have up to a billion elements (with the 
limiting factor from running larger runs being memory; all the simulations reported here use 
< 64 Gb). Our typical simulation involves 10® sheets with 10® PM elements, which is more 
than sufficient for the simulation to estimate the nonlinear P{k) to 1% accuracy. In fact, the 
smaller number of resolution elements needed is a significant advantage of ID studies: In ID, 
a simulation volume of 


Did 


1.3 X 10^ Mpc 


^_2 f 0.01 Mpc 

h% 


(A.l) 


is required to be able to measure the power spectrum with accuracy of fi% (where fi% is in 
units of 1%), assuming Gaussian statistics and a bin width in A;-space of Ak. In contrast, in 
3D a volume of A: = |k| requires 
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4 X 10® Mpc® 


f-2 ( 0-1 Mpc 

I k 


-1 


0.01 Mpc 

aF 
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(A.2) 


assuming all wavenumbers that fall within a shell of width Ak are binned to measure the 
power spectrum. Even if convergence is reached with merely 1 particle per Mpc""- a typical 
resolution of 3D simulations aimed at precisely measuring the evolution at mildly nonlinear 
scales [74] and what we find here is needed in ID -,4x 10® particles are required in dimension 
n = 3 for k = 0.1 Mpc“® and Ak = 0.01 Mpc“® versus just 1 x lO'^ in n = 1 at any k and the 
same Ak to reach 1% fractional precision. 

Figure 8 shows convergence tests for the z = 0 GDM-like case. The lefthand panels 
show the convergence in PM grid size (denoted “PM”), number of mass elements (“N”), box 
size (“Box”), where each estimate of P{k) has used 30,000 modes. Runs with the same box 
size have the same random numbers for overlapping modes. We find that the most important 
variable is resolution, with convergence at A: < 1 Mpc“® at the percent level once the PM 
resolution is < 1 Mpc and the mean inter-sheet spacing satisfies the same inequality. In the 
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main body of the paper, we use a fiducial resolution of 0.1 Mpc, box sizes > 10^ Mpc, and 
number of particles and PM grids > 10® Mpc. The inequalities arise because in the CDM- 
like case our investigations beneht from even higher precision than ~ 1% and so we run the 
largest simulation possible (for the power-law cases equality holds). The righthand panel in 
figure 8 compares the power spectrum in simulations that change the initialization redshift 
and time stepping. For the fiducial choices of z = 20 and Ay^ = 0.001 the calculations are 
converged to a few tenths of a percent level in both parameters. (That the lower initialization 
redshift simulation has similar residuals to the small time step simulation suggests that much 
of the error occurs at the highest redshifts, likely indicating that our time stepping criterion 
becomes too coarse there.) We conclude that our code is converged to sub-percent for the 
desired calculations in the CDM-like cosmology. 

We have run separate convergence tests for the power-law cosmologies discussed in 
section 3 that vary the initialization redshift (but keep 0.1 Mpc resolution and the fiducial 
time step). For these calculations we maintain the box sizes and resolutions indicated above. 
Only for the n = 2 power-law case do we have to increase the initialization redshift to z = 40 
in order to achieve the desired level of convergence. 


B The relationship between SPT and LPT 

B.l Proof that ID LPT is identical to ID SPT at every order in <5^ 

We first derive a useful recurrence relation obeyed by both SPT and LPT. Starting with the 
SPT relations (copying eq. 2.10): 

Fn = (2n -|- l)Xn Yn] Gn = 3An -|- nYrii (B.l) 

and solving for Y^ using the first equation and, then, substituting this relation into the second 
equation yields 


— 3A„ -|- n [Fn — (2n -|- l)Xn\ , 

(B.2) 

= nFn - (2n + 3){n- 1)X„, 

(B.3) 

n—1 , 

— n / j n—rri’i 

1 -^1 
m=l 

(B.4) 


where the last line uses the definition of Xn (eq. 2.11). 

We now show that, after symmetrizing, the SPT relation given by eq. (B.4) is identical 
to the symmetrized LPT recursion relation. Our starting point is to derive an equation 
that relates velocities and overdensities in LPT (which turns out to be the same recurrence 
relation as just derived, eq. B.4). This relation will allow us to derive Gn from the Fn of the 
Zeldovich approximation and then to show that the two satisfy both SPT recursion relations. 
We can relate the two using the following equation for the LPT momentum field: 

(1 Sza)uza = y -Q- ^ ZA ( g ')]- (B.5) 


Going to Fourier space and taking the gradient of the above yields the relation 
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Identifying terms in this equation that have the same power in 6l yields 

^ZA*f^ZA ^^ZA^'^A 

/■ s /■ “S 

Gn ^ ^ GrriFn—m “t“ ^ ^ GmFn—m — (B-7) 

m=l m=l 

where we have labeled the less-trivial convolution terms, and are using the shorthand Gn = 
Gn{ki, ■■■ ,kn) and Fn-m = Fn{kn, • • • , km)- Simplifying eq. (B.7) gives 

Gn + -f^GmFn-m = nFn, (B.8) 

1 Ai 
m=l 

and we can easily symmetrize this relation by averaging each term by taking all permutations 
of the ki- We will denote this symmetrizing operation on term X as = Sym[X]. Note 

that the symmetrized version of eq. (B.7) is the same relation as satisfied by SPT (eq. B.4). 
Even though in both LPT and SPT we derived this recursion relation using our ID expres¬ 
sions, this derivation generalizes to 3D and shows that SPT and LPT share at least one 
recursion relation. Physically this relation owes to the continuity of u and 6. 


We now show that = F^'^ satisfies recursion relation (B.8), using our result that 
^sym _ ^ 2.34). For this to hold, eq. (B.8) reduces to 
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(B.9) 

(B.IO) 

(B.ll) 

(B.12) 

(B.13) 

(B.14) 

(B.15) 


where to reach eq. (B.13) we have used that Sym[7Li] = {m/n)k and Sym[iL 2 ] = {[n — m\/n)k. 
The final line in this reduction (eq. B.14) is an identity (which we have verified up to n = 1000 
using Mathematica and which holds in the 1 Sterling’s approximation limit). 

Now that we have F^'^ and Gn^™ in ID LPT, it is simple to show that they also satisfy 
the SPT recursion relations. Since 
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noting that k?/{K 1 K 2 ) = k/Ki + k/K 2 , it follows that 

(B.17) 

" 2n + 3"’" 271 + 3^^’ ^ ^ 

where Xn and Yn were defined in eqs. (2.11) and (2.12). It is trivially seen that these satisfy 
the symmetrized SPT recursion relations, eq. (2.10). ■ 

B.2 Derivation of SPT from LPT 

To generate SPT-like expressions for P{k) from LPT, Plpt(^) can be expanded in powers of 
PL{k): 

PhPTik) = j 0 -^( 9 ) + -) > (^-18) 

= PL + \j dqP’^^Vla\q) + 0{Pl), (B.19) 

6 {[aX? + 8 (k^]'k^]'") + 2 + O(P|),(B.20) 
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(B.21) 


where the second line used that which integrating by parts equals e*^'?Vq 

up to surface terms, and the last line substitutes for cj^ using eq. (2.35). As anticipated, the 
last line in the above equation is identical to what we found for 1-loop SPT (eq. 2.13), and 
Appendix B.l proves that this correspondence holds at all loop orders. 

Similarly, we can also start with the correlation function in the LPT (eq. 2.38) and 
expand: 
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where (• • • )o- denotes an average over x' with a Gaussian kernel of width a. The last line 
allows us to evaluate each term in the expansion using Wick’s theorem = (m— 1 )!!ct^™', 

which reduces to 

(B.23) 

\x / 

,2 

(B.24) 


izA{x) = ^ + + (P4['^"l"') + -y\[<AT + O (^) 

j-ao 2 

= iL{x) + S^Lixf + 4^l(x)' j dxiiix) + —iL{x)'' + 


o[^- 


Primes denote derivatives with respect to the argument, and for the error term we have 
treated, e.g., ~ equivalent. In addition, we have used that /‘^■ 

Eq. (B.24) is identical to our 1-loop expression for the correlation function (eq. 4.1), and it 
is the Fourier transform of eq. (B.21). 
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B.3 Stochastic fluctuations in standard theories 

At large scales it is well known that the nonlinear evolution of a power spectrum that has no 
power as fc —)• 0 is to develop a tail in which P{k) ~ plus terms that are higher in power 
[75]. This limiting behavior can also be derived from our SPT and LPT expressions. We 
have shown that the power spectrum can be expanded to yield (e.g., eq. B.18) 

Pi-ioop(fc) = y^-yCj2(g) + y((Tf+ o-^)2(g) + --y , (B.25) 

= Piik) + y (cjf + cj|) * {af + af){k) + • • • , (B.26) 

where we have written explicitly the terms up to one-loop order. The purely stochastic term, 
oc k^a'^ * ag , scales as as A; —?• 0 as this term is convolving two broadband high-pass filtered 
fields (e.g., the convolution does not depend on k in the limit k <^k'). This behavior is not 
as apparent in P 22 in the text (i.e. eq. 2.15 or equivalently eq. B.21). 


C Equations in the effective field theory of large-scale structure 
C.l Derivation of smoothed equations of motion 


The presentation in this section follows the derivation of effective equations of motion in 
[21] and [31], but also comments on some recent improvements to these original derivations. 
We start with the collisionless Boltzmann equation for the particle distribution function 
f{p, r), where p is the momentum including the Hubble flow contribution and r is the proper 
coordinate: 


dfir,p,t) 

dt 




mVr^Vpf = 0 , 


(C.l) 


where includes the homogeneous component so, e.g., /A in Einstein de-Sitter. We 

smooth this equation on a comoving scale A“^, again using the shorthand Xi = f dr'Wf^{r — 
r')X{r') for the smoothed long wavelength field. Since our smoothing is a linear operation, 


^ + -^rfi =m [ dPW^ir - P)Xr'Hr')Xpf{r',p). 
at m J 

We define the density, momentum, and velocity dispersion fields to be 


(C.2) 


P = mjdpf{r,p,t), (C.3) 

TT = j dppf{r,p,t), (C.4) 

j dpp‘^f{r,p,t). (C.5) 


To proceed we take moments of equation (C.2) by integrating each term over f dpp^. Since 
the right hand side of equation (C.2) is a total derivative with respect to p, the zeroth moment 
yields the continuity equation: 


dtPi + '^rT^i = j dr'W{r — r')[Vr'4>(r')] / dpVpf{r',p) = 0, 


(C.6) 
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where in the last step we have used that the integral of a derivative is zero. If we define 
v’l = TTi/Pi (the superscript on vf is to denote that this is not truly a long-wavelength field), 
we obtain the continuity equation in its usual form: 


Pl+Vr{piv'^)=0. (C.7) 

Next consider the 1®* moment, which should give an Euler-like or Navier-Stokes-like 
equation. Multiplying by p and integrating Eq. (C.2) over p gives 

dtiTi + Vra‘^ i = m j dr'W{r - r')[Vr'^{r')] J dppdpf{r',p), (C.8) 

where ci^ ^ is the smoothing of , not the square of the smoothed sigma. Integrating by 
parts the integral over dp yields 

dtni + Vrcrh = - j dr' IE(r - r')[V,/$(r')]p(r'). (C.9) 

Writing vr; = piv^ and, then, pulling out the linear terms from the RHS, we have 

pidtvf + v\dtpi + VrCF^ i = - j dr' W{r - r')[Vr'^{r')]p{r'), (C.IO) 


-Vr^ipi + pVAA -p I dr' W{r - r')[Vr'(t){r')]5{r'). (C.Il) 

To reach this equation, we assumed that the homogeneous part of <h, d>/j(r), can be ap¬ 
proximated as being evaluated at r rather than r' in order to pull it out of the inte¬ 
gral. To calculate the error from this approximation, we expand around r' as = 

+ Vrd>/j(r)(r^ — r) + 2 V^<h/i(r)(r' — r)^ -|- .... since it is smoothly varying (in Einstein 
de-Sitter and ACDM this expansion is exact at second order). Eurthermore, we specialized to 
the Gaussian functional from for IE(r — r') = A exp[—A^(r — r')^/2]/\/^ and to an Einstein 
deSitter cosmology such that 

p j dr' W{r - r') [V^,$,,(/) - V^$/j(r)] 6{r') = 


where A is comoving. Thus, the error from approximation is suppressed by A^, justifying our 
approximation. 

Starting again with eq. (C.ll), we can use the continuity equation to eliminate dtpi, 
which yields: 


p J dr' W{r — r')[V‘^^h{r)ir'— r)]5{r'), 
—AnGp^ f dr' W{r — r')[r' — r]5{r'), 


-H^p / dr' A-^VrW{r - r')6{r'), 
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,yAr') 
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(C.I2) 


dtV^ - Pi ^v\Vr{Plv\) + Pi = -Pi ^ {p [VrCpdji - pV■ (C.I3) 
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Furthermore, we can use the Poisson equation, = —ATrGpS, to eliminate 6 on the RHS: 

dtvf + v\Vrv\ - pY^Vr{pivf) + pJ^V ^ - V , 
or 


dtv\ + v\\/ rv\ + + 



1 

SttG/Oz 


Vr ([Vr<^Vr<^]; 


V r4^l^r4^l) ) 


which can be written as 

dtv'l + v\Vrv\ + Vr^l = --V^TA. (C.14) 

Pi 

We decompose the long-wavelength stress tensor as t\ = ka + Sa, the sum of a kinematic 
term, 

- Pivf, (C.15) 

and a gravitational term. 


Sa = {[Vr(pVr(f)]l - Vr(plVr4>l) ■ (C.16) 

The above equations are the dynamical equations of EFTLSS. To make better contact 
with the velocity that appears in standard perturbation theories as well as to be able to 
measure the coefficients of the effective theory on a grid, we formulate the EFTLSS in terms 
of the Eulerian velocity vi (rather than uj*) via an expansion in dp. 

v\^- = = vi + - (Hz - vim + 0{6f), (C.17) 

Pi l + di 

= vi + [u(,5 - <50]z(l - Si) + Oi5f,A-^). (C.18) 


We note that the Eulerian velocity vi is still sensitive to short wavelengths as detailed in [22] , 
affecting > 1 loop order predictions (which are not of concern here). 

Using this expansion, the continuity equation becomes, to linear order in the long- 
wavelength fields. 

Pi + ^ripm) = -Vr {pi[v{6 - Si)]i - p{[v6]a - viSi)Si} + 0{6f), (C.19) 

= -p v,[u(5 - + 0{6f, A-2). (C.20) 


Where the latter relation holds to the extent that large-scale fields do not affect small-scale 
averaging (i.e. it drops terms suppressed by A“^). Nicely, we now have terms on the RHS 
of the continuity equation that in effect generate terms that [23] showed are not forbidden 
by symmetry and that are needed to cancel divergences that appear in the 1-loop velocity 
power spectrum or higher loop calculations. 

Turning to the Euler equation and plugging in our expansion for u]': 

dtvi+viVrVi+S/r^i = -p~^VrTA-dt{[v{6-5i)]i{l-6i)}-S/r {vi[v{6 - 6i)]i{l - 6i))+0{Sf,A~'^). 

(C.21) 

Note that, e.g., [Vu(5 — (^/)]a can contribute a constant and so we maintain terms that are 
higher order. The above physical coordinate equations can all be written as the comoving 
coordinate equations that are solved perturbatively in SPT (e.g., [75]; Vr —)• a“^V and 
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dt = dt — Hx'V). The only difference with respect to SPT is that the terms with and 
— Si)]i. These additional terms are sensitive to short-wavelength perturbations that are 
nonlinear in the cosmologies we are considering. 

What we need to solve for the nonlinear density is just the nonlinear equation for 6, 
which plugging in our new Euler and continuity equations and going to comoving coordinates 
yields 

-a^'H^{a)dl5i - a {2%^{a) + a'H{a)d'H{a)/da) da5i + ^TrGpa^di = {aUda + n) V{diui) - V{uiVui) 

-P-^VVa + Vanda{[u{6 - 6i)]i6i) - ([n(5 - 6i)]iui) - Vn[ui6 - <5z)];(l - <5;) +0{6f,A-^). 

^-V-^ 

terms not in SPT 

To solve the EFTLSS equations perturbatively, the new short-wavelength sensitive terms 
are expanded in terms of the long-wavelength fields (at all previous times; [22, 24]). The coef¬ 
ficients of these terms are then estimated from a simulation in Appendix D. Another approach 
to deriving the effective equations is to include all possible terms allowed by the symmetries of 
the problem - homogeneity &: isotropy, Galilean invariance, mass and momentum conserva¬ 
tion, and the equivalence principle - as “source” terms to the continuity and Euler equations 
[23]. A quick summary in the spirit of this approach is described in the ensuing appendix. 

C.2 Symmetries approach 

In section 2.4 we described one method for overcoming the deficiencies of SPT. In that 
approach, an attempt was made to construct dynamical equations for the smoothed matter 
overdensity and velocity fields. An alternative approach [23], which arrives at the same extra 
terms as we find in eq. (2.51) is to write down all of the terms consistent with the symmetries. 

(We note that the symmetries are also taken account in the first method; this method is just 
more pedagogical.) We briefly review this approach here. 

When computing ensemble averages of cosmological fields we can imagine taking the 
average in two steps. First we average over the small-scale (and potentially nonlinear) fluctu¬ 
ations and then over the longer wavelength modes (for which we assume perturbation theory 
is applicable). In the intermediate stage, after “integrating out” the small-scale perturba¬ 
tions, we have a set of evolution equations for the long-wavelength modes in which the effects 
of the small-scale perturbations are encoded as additional terms beyond the terms one usually 
encounters in perturbation theory. If we write the effective equations for the long-wavelength 
perturbations with the (standard) linear terms of the left-hand side (e.g. eqs. 2.41 and 2.42), 
then the right hand side must contain all other terms allowed by the symmetries of the prob¬ 
lem, though with unknown coefficients. In addition to the standard quadratic terms [V((5u) 
and V(mVu)] we must handle extra terms such as those in ta of eq. (2.42). One can write 
all such terms consistent with the symmetries to a given order in 6, u, and derivatives, with 
the full effective equations having infinitely many terms. 

To make progress these equations can be solved perturbatively, expanding in the linear 
overdensity, 6l, and only treating terms in the equations at or below the desired order 
in 5, u, and derivatives. We arrange the source terms in a derivative expansion (i.e. in 
Fourier space an expansion in positive powers of k'^) and also a Taylor expansion in powers 
of the long-wavelength perturbations (upon which the small-scale dynamics can depend). 

At lowest order, the only dependence can be on 5l since at linear order Ol = 5l and the 
vorticity vanishes. At the next nontrivial order and in Fourier space the allowed terms go as 
a constant times - it cannot be a constant times 6i or mass would not be conserved nor 
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Figure 9. Contribution of stochastic (small-scale) fluctuations to the power spectrum at z = 1 for our 
CDM-like case. The stochastic term is computed by differencing 6 in two simulations that have the 
same modes with k < 2 ttA and different modes at larger fc, and then dividing by two. Interestingly, 
the curves quickly reach the fc® small wavenumber limit once k < 27rA. The dashed lines show the fc® 
asymptote that these curves should reach as fc —> 0. The stochastic term is controversially dropped in 
previous EFTLSS calculations, but the quick asymptote of this term to ^ k^ suggests that dropping 
it is indeed an excellent approximation. The dot dashed curve shows a possible upper bound on Pj 
discussed in the text. 


times k6 l or momentum conservation would be violated - or are independent of 5 l with RMS 
scaling as /c^ due to mass and momentum conservation [75]. These terms will have unknown 
coefficients with unknown time-dependence. In the main body, we justify why these terms 
should be thought of as third order terms. When integrated against the Green’s function 
these new terms generate contributions to the density going as (5/-independent (J) or as a 
coefficient times k‘^ times the linear theory overdensity. 

When computing the power spectrum of 5i the first term only comes in as Pj in 
eq. (2.51). We know that Pj goes as as /c —?■ 0 but the behavior outside this limit is 
unknown. The other term contributes k‘^PL{k) to lowest order, with an unknown coefficient 
(2 q:c,a in eq. 2.51) that is typically fit for. The final result for the density power spectrum is 
the ‘usual’ perturbative solution plus a stochastic term (going as k^ as k ^ 0) plus a term 
going as k‘^PL{k). This is identical to eq. (2.51). 

D Importance of stochastic terms 

Previous calculations of P{k) in EFTLSS have dropped the stochastic term Pj in equa¬ 
tion (2.51) - the term that does not correlate with any term in the expansion in the smoothed 
density and velocity - because of its weaker fc —>■ 0 scaling [22, 31]. This is curious, as halo 
models for large-scale structure explain nonlinear evolution with just a stochastic term (scal¬ 
ing with the abundance of halos) that goes as k^ [33-35]. That halo models are at least 
qualitatively successful suggests that the stochastic term may have a non-trivial dependence 
on wavenumber at mildly nonlinear wavenumbers, which may invalidate our argument for its 
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smallness based purely on the /c —)• 0 limit. This objection was raised in [37, 38], which used 
the halo profile to motivate a completely different expansion from that in EFTLSS for addi¬ 
tional contributions to P{k) (although we note that their corrections may not be “stochastic” 
in our sense, i.e. meaning that they do not correlate with the smoothed fields). Naively, the 
impact of the stochastic term and the profiles of nonlinear systems on P{k) should be at least 
as important in ID as in 3D because (1) even though their profile may be different, virialized 
systems still should have the same extent (half the turnaround radius),and (2) there are 
fewer realizations of the small-scale modes in ID than in 3D. 


However, we find dropping the stochastic term to be an excellent approximation (and 
accurate at the percent level a factor of ~ 3 beyond the smoothing scale). Figure 9 shows 
the contribution of stochastic (small-scale) fluctuations that are ignored in FFTLSS to the 
z = 1 matter power spectrum for our CDM-like case. The stochastic power, Pj, is computed 
by differencing the 6 in simulations that have the same random numbers for modes with 
k < 2 ttA and different ones for k > 27rA, and then dividing by two. We note that A in this 
case is not the smoothing scale but rather the coherence scale of the two simulations. The 
quantity that enters into the FFTLSS power spectrum for finite A, Pj^\ is essentially W\ 
times what is plotted at /c < 27rA. Beyond 27rA (the filled circles), there is an intermediate 
range of scales where the stochastic terms scale as k in kP{k) that is most apparent for the 
smaller A values in figure 9 (as in the halo model). However, the curves quickly reach the k^ 
large-scale limit, and hence the halo model ansatz that the shape of halos are an important 
contribution to the mildly nonlinear power is not correct, at least in ID. 


The Pj curves in figure 9 do not give an estimate for the actual value for Pj^\ as A —)> oo. 
Instead, it is how quickly these curves reach the k^ asymptote that is of most interest, as 
the steep asymptotic scaling is what is used to drop these terms in current applications 
of EFTLSS. In EFTLSS these stochastic contributions should sum, at lowest order, with 
UV sensitive terms in P 22 ,A to generate a A-independent term, but their functional form 
potentially matters near A. A way to place an upper bound on Pj as A —)• oo (to the extent 
that Pj in this limit makes sense at all) is to take the square of the difference between 
the nonlinear density field and the Zeldovich approximation density field, normalizing by 
a transfer function so that Zeldovich has been rescaled to have the same power spectrum: 
i.e., to calculate {{5nl — T5za)‘^) where T = {Pnl/PzaY^'^ and 5ml and Pml are the 
nonlinear density and power spectrum computed in the simulation. This quantity, which also 
equals PnlY ~ where r is the cross correlation coefficient between 5ml and 5zAi is the 
part of the nonlinear density that does not correlate with the Zeldovich density and hence 
more or less excludes the SPT terms at all orders. It is comprised of both the stochastic 
term but also non-stochastic contributions from effective terms that do not trace 5za- This 
quantity is shown by the dot-dashed curve, which is down by three orders of magnitude from 
the simulation power at A: = O.I Mpc“^. The same statistic is similarly suppressed in 3D 
[53]. 
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Figure 10. Estimates for 2ac,A constructed by first estimating from nonlinear simulations (using 
eq. D.l) and, then, assuming that oc a to solve for 2ac,A via eq. (2.50). The solid blue curve 
corresponds to the 2ac,A needed to fit the power spectrum in our simulations, using the 2ac estimated 
in section 2.4.3. The points with error bars show the estimated values from Xj^. All curves and points 
are computed for z = 0.5. 


D.l Estimating the new coefficients of the effective theory 

Equation (2.43) for Xa provides the dependence of the smoothed fields in terms of a “stress 
tensor” that is sensitive to small scales, which we then expanded in terms of large-scale modes 
and effective coefficients. By measuring X\ and how it correlates with large-scale modes, we 
can measure the coefficients of this theory and thus not have to fit free parameters when we 
compare the EFTLSS predictions to the power spectrum in simulations. An estimator for 
the relevant linear combination of these coefficients, from the small-scale behavior of a 
simulation is ^ 

p^t= Y.^kXA{k)6li{k), (D.l) 

\ k / k 

where the sum should be over fe <C A and we choose the weightings, Wk, that maximize the 
S/N assuming Gaussianity. We compute X\(k) from the Eulerian velocity and overdensity 
field on our ID simulation grid. The new EFTLSS parameter that enters into the 1-loop 
power spectrum, acA, is then related to Ctot,A via eq. (2.49). 

This procedure has been found to produce the 2ac needed in [31] using a somewhat 
different method to the above estimator, although they did choose what seemed like a very 
fine-tuned range of separations of 7—10 Mpc for A“^ = 3 Mpc (13—18 Mpc for A“^ = 6 Mpc). 
A variant of this procedure where the time evolution of is measured explicitly (without 
assuming a long-wavelength expansion) was also found to work in [62]. 

^^Virialized systems are much less dense in ID than halos in 3D. However, we do not expect that the density 
in the virialized region matters for the power spectrum at mildly nonlinear scales. Imagine taking a kernel 
and convolving the density field on a scale of a few Mpc (small enough to not affect mildly nonlinear scales). 
This dramatically reduces the density contrast in the field. However, the power spectrum on mildly nonlinear 
scales is not affected since it is multiplied by the Fourier transform of the smoothing kernel squared. Thus, 
for the matter evolution on mildly nonlinear scales it does not matter what the structure of systems is on less 
than a few Mpc scales, larger than the sizes of all halos in 3D concordance cosmology. 
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Figure 10 shows our attempts to measure from a ID simulation in our CDM-like 
cosmology at z = 0.5 for different A, where what is plotted is 2ac,A, which is proportional to 
The proportionality constant is determined using eq. (D.l) and oc a as required to 
scale with Pl3, but it is likely the true temporal scaling is somewhat different and so these 
estimates are only approximate. The solid curve corresponds to the 2 q:c,a required to match 
the power spectrum in simulations, determined in section 2.4.3. The points with error bars 
show the estimated value from Xa (red crosses). These points are in agreement with the 
simulation fit. The error bars are rather large, which owe to difficulties in estimating the 
Eulerian space velocity that appears. (Indeed, these calculations required a higher density 
of sheets per cell than our fiducial simulations for the error bars do be constraining at all. 
This computation was done with 10® sheets and lO'^PM grid cells in a 10® Mpc box.) It likely 
reduces the variance to formulate Xa in terms of the momentum, as done in [31]. 

References 

[1] Y. B. Zel’dovich, Gravitational instability: An approximate theory for large density 
perturbations.^ A&A 5 (Mar., 1970) 84-89. 

[2] M. Davis, G. Efstathiou, C. S. Frenk, and S. D. M. White, The evolution of large-scale 
structure in a universe dominated by cold dark matter, ApJ 292 (May, 1985) 371-394. 

[3] V. Springel, S. D. M. White, A. Jenkins, C. S. Frenk, N. Yoshida, L. Gao, J. Navarro, 

R. Thacker, D. Groton, J. Helly, J. A. Peacock, S. Gole, P. Thomas, H. Couchman, A. Evrard, 
J. Golberg, and F. Pearce, Simulations of the formation, evolution and clustering of galaxies 
and quasars, Nature 435 (June, 2005) 629-636, [astro-ph/0504097]. 

[4] M. Tegmark et al.. The Three-Dimensional Power Spectrum of Galaxies from the Sloan Digital 
Sky Survey, ApJ 606 (May, 2004) 702-740, [astro-ph/0310725]. 

[5] S. Cole et al., The 2dF Galaxy Redshift Survey: power-spectrum analysis of the final data set 
and cosmological implications, MNRAS 362 (Sept., 2005) 505-534, [astro-ph./0501174]. 

[6] W. Hu and N. Sugiyama, Small-Scale Gosmological Perturbations: an Analytic Approach, ApJ 
471 (Nov., 1996) 542, [astro-ph/9510117]. 

[7] Planck Collaboration, P. A. R. Ade, N. Aghanim, C. Armitage-Caplan, M. Arnaud, 

M. Ashdown, F. Atrio-Barandela, J. Aumont, C. Baccigalupi, A. J. Banday, and et al., Planck 
2013 results. XVI. Gosmological parameters, A&A 571 (Nov., 2014) A16, [arXiv; 1303.5076]. 

[8] A. G. Sanchez et al.. The clustering of galaxies in the SDSS-III Baryon Oscillation 
Spectroscopic Survey: cosmological implications of the large-scale two-point correlation 
function, MNRAS 4:26 (Sept., 2012) 415-437, [arXiv; 1203.6616]. 

[9] M. H. Goroff, B. Grinstein, S.-J. Rey, and M. B. Wise, Coupling of modes of cosmological mass 
density fluctuations, ApJ 311 (Dec., 1986) 6-14. 

[10] N. Makino, M. Sasaki, and Y. Suto, Analytic approach to the perturbative expansion of 
nonlinear gravitational fluctuations in cosmological density and velocity fields, PRD 46 (July, 
1992) 585-602. 

[11] B. Jain and E. Bertschinger, Second-order power spectrum and nonlinear evolution at high 
redshift, ApJ 431 (Aug., 1994) 495-505, [astro-pli/9311070]. 

[12] T. Buchert, Lagrangian theory of gravitational instability of Friedman-Lemaitre cosmologies 
and the ’Zel’dovich approximation’, MNRAS 254 (Feb., 1992) 729-737. 

[13] T. Buchert, Lagrangian Theory of Gravitational Instability of Friedman-Lemaitre Cosmologies - 
a Generic Third-Order Model for Nonlinear Clustering, MNRAS 267 (Apr., 1994) 811, 
[astro-ph/9309055]. 


- 40 - 


[14] P. Catelan, Lagrangian dynamics in non-fiat universes and non-linear gravitational evolution, 
MNRAS 276 (Sept., 1995) 115-124, [astro-pli/9406016]. 

[15] J. Ehlers and T. Buchert, Newtonian Cosmology in Lagrangian Formulation: Foundations and 
Perturbation Theory, General Relativity and Gravitation 29 (June, 1997) 733-764, 
[astro-ph/9609036]. 

[16] F. Bernardeau, S. Colombi, E. Gaztanaga, and R. Scoccimarro, Large-scale structure of the 
Universe and cosmological perturbation theory, physrep 367 (Sept., 2002) 1-248, 

[astro-ph/0112551]. 

[17] P. McDonald, Dark matter clustering: A simple renormalization group approach, PRD 75 
(Feb., 2007) 043514, [astro-pli/0606028]. 

[18] M. Crocce and R. Scoccimarro, Renormalized cosmological perturbation theory, PRD 73 (Mar., 
2006) 063519, [astro-ph/0509418]. 

[19] S. Matarrese and M. Pietroni, Resumming cosmic perturbations, JCAP 6 (June, 2007) 26, 
[astro-ph/0703563] . 

[20] T. Matsubara, Resumming cosmological perturbations via the Lagrangian pieture: One-loop 
results in real space and in redshift space, PRD 77 (Mar., 2008) 063530, [arXiv:0711.2521]. 

[21] D. Baumann, A. Nicolis, L. Senatore, and M. Zaldarriaga, Cosmological non-linearities as an 
effective fluid, JCAP 7 (July, 2012) 51, [arXiv: 1004.2488]. 

[22] J. J. M. Carrasco, S. Foreman, D. Green, and L. Senatore, The Effective Field Theory of Large 
Scale Structures at two loops, JCAP 7 (July, 2014) 57, [arXiv: 1310.0464]. 

[23] L. Mercolli and E. Pajer, On the velocity in the Effective Eield Theory of Large Scale 
Structures, JCAP 6 (Mar., 2014) 6, [arXiv: 1307.3220]. 

[24] S. M. Carroll, S. Leichenauer, and J. Pollack, Consistent effective theory of long-wavelength 
cosmological perturbations, PRD 90 (July, 2014) 023518, [arXiv: 1310.2920]. 

[25] R. A. Porto, L. Senatore, and M. Zaldarriaga, The Lagrangian-space Effective Eield Theory of 
large scale structures, JCAP 5 (May, 2014) 22, [arXiv: 1311.2168]. 

[26] N. Padmanabhan and M. White, Calibrating the baryon oscillation ruler for matter and halos, 
PRD 80 (Sept., 2009) 063508, [arXiv:0906.1198]. 

[27] J. Carlson, B. Reid, and M. White, Convolution Lagrangian perturbation theory for biased 
tracers, MNRAS 429 (Feb., 2013) 1674-1685, [arXiv: 1209.0780]. 

[28] M. White, The Zel’dovich approximation, MNRAS 439 (Apr., 2014) 3630-3640, 

[arXiv: 1401.5466]. 

[29] D. J. Eisenstein, H.-J. Seo, and M. White, On the Robustness of the Acoustic Scale in the 
Low-Redshift Clustering of Matter, ApJ 664 (Aug., 2007) 660-674, [astro-ph/0604361]. 

[30] S. Tassev, Lagrangian or Eulerian; real or Fourier? Not all approaches to large-scale structure 
are created equal, JCAP 6 (June, 2014) 8, [arXiv: 1311.4884]. 

[31] J. J. M. Carrasco, M. P. Hertzberg, and L. Senatore, The effective field theory of cosmological 
large scale structures. Journal of High Energy Physics 9 (Sept., 2012) 82, [arXiv: 1206.2926]. 

[32] L. Senatore and M. Zaldarriaga, The IR-resummed Effective Field Theory of Large Scale 
Structures, ArXiv:IfOf-SJSf (Apr., 2014) [arXiv: 1404.5954]. 

[33] U. Seljak, Analytic model for galaxy and dark matter clustering, MNRAS 318 (Oct., 2000) 
203-213, [astro-ph/0001493]. 

[34] C.-P. Ma and J. N. Fry, Deriving the Nonlinear Cosmological Power Spectrum and Bispectrum 
from Analytic Dark Matter Halo Profiles and Mass Functions, ApJ 543 (Nov., 2000) 503-513, 
[astro-ph/0003343] . 


- 41 - 


[35] A. Cooray and R. Sheth, Halo models of large scale structure, physrep 372 (Dec., 2002) 1-129, 
[astro-ph/0206508]. 

[36] N. Afshordi, How well can (renormalized) perturbation theory predict dark matter clustering 
properties?, PRD 75 (Jan., 2007) 021302, [astro-ph/0610336]. 

[37] I. Mohammed and U. Seljak, Analytic model for the matter power spectrum, its covariance 
matrix and baryonic effects, MNRAS 445 (Dec., 2014) 3382-3400, [arXiv: 1407.0060]. 

[38] U. Seljak and Z. Vlah, Halo Zeldovich model and perturbation theory: dark matter power 
spectrum and correlation function, ArXiv:1501.07512 (Jan., 2015) [arXiv: 1501.0751]. 

[39] C. D. Rimes and A. J. S. Hamilton, Information content of the non-linear power spectrum: the 
effect of beat-coupling to large scales, MNRAS 371 (Sept., 2006) 1205-1215, 

[astro-ph/0511418]. 

[40] B. D. Sherwin and M. Zaldarriaga, Shift of the baryon acoustic oscillation scale: A simple 
physical picture, PRD 85 (May, 2012) 103523, [arXiv: 1202.3998]. 

[41] T. Matsubara, The Correlation Function in Redshift Space: General Formula with Wide-Angle 
Effects and Cosmological Distortions, ApJ 535 (May, 2000) 1-23, [astro-ph/9908056]. 

[42] N. E. Chisari and M. Zaldarriaga, Connection between Newtonian simulations and general 
relativity, PRD 83 (June, 2011) 123505, [arXiv: 1101.3555]. 

[43] E. T. Vishniac, Why weakly non-linear effects are small in a zero-pressure cosmology, MNRAS 
203 (Apr., 1983) 345-349. 

[44] R. Scoccimarro and J. Frieman, Loop Corrections in Nonlinear Cosmological Perturbation 
Theory, ApJS 105 (July, 1996) 37, [astro-ph/9509047]. 

[45] J. J. M. Carrasco, S. Foreman, D. Green, and L. Senatore, The 2-loop matter power spectrum 
and the IR-safe integrand, JCAP 7 (July, 2014) 56, [arXiv: 1304.4946]. 

[46] A. J. S. Hamilton, C. D. Rimes, and R. Scoccimarro, On measuring the covariance matrix of 
the non-linear power spectrum from simulations, MNRAS 371 (Sept., 2006) 1188-1204, 
[astro-ph/0511416]. 

[47] M. Takada and W. Hu, Power spectrum super-sample covariance, PRD 87 (June, 2013) 

123504, [arXiv: 1302.6994]. 

[48] A. Meiksin and M. White, The growth of correlations in the matter power spectrum, MNRAS 
308 (Oct., 1999) 1179-1184, [astro-ph/9812129]. 

[49] R. de Putter, C. Wagner, O. Mena, L. Verde, and W. J. Percival, Thinking outside the box: 
effects of modes larger than the survey on matter power spectrum covariance, JCAP 4 (Apr., 
2012) 19, [arXiv: 1111.6596]. 

[50] T. Baldauf, U. Seljak, L. Senatore, and M. Zaldarriaga, Galaxy bias and non-linear structure 
formation in general relativity, JCAP 10 (Oct., 2011) 31, [arXiv: 1106.5507]. 

[51] N. Padmanabhan, M. White, and J. D. Cohn, Reconstructing baryon oscillations: A 
Lagrangian theory perspective, PRD 79 (Mar., 2009) 063523, [arXiv:0812.2905]. 

[52] Y. Noh, M. White, and N. Padmanabhan, Reconstructing baryon oscillations, PRD 80 (Dec., 
2009) 123501, [arXiv:0909.1802]. 

[53] S. Tassev and M. Zaldarriaga, The mildly non-linear regime of structure formation, JCAP 4 
(Apr., 2012) 13, [arXiv: 1109.4939]. 

[54] S. F. Shandarin and Y. B. Zeldovich, The large-scale structure of the universe: Turbulence, 
intermittency, structures in a self-gravitating medium. Reviews of Modern Physics 61 (Apr., 
1989) 185-220. 


- 42 - 


[55] V. Sahni and P. Coles, Approximation methods for non-linear gravitational clustering, physrep 
262 (Nov., 1995) 1-135, [astro-ph/9505005]. 

[56] B. Grinstein and M. B. Wise, On the validity of the Zel’dovich approximation, ApJ 320 (Sept., 
1987) 448-453. 

[57] C. Rampf and T. Buchert, Lagrangian perturbations and the matter bispectrum I: fourth-order 
model for non-linear clustering, JCAP 6 (June, 2012) 21, [arXiv; 1203.4260]. 

[58] A. N. Taylor, Statistics of self-gravitating fluctuations, in Cosmic Velocity Fields (F. Bouchet 
and M. Lachieze-Rey, eds.), p. 585, 1993. 

[59] P. Schneider and M. Bartelmann, The power spectrum of density fluctuations in the Zel’dovich 
approximation, MNRAS 273 (Mar., 1995) 475-483. 

[60] K. B. Fisher and A. Nusser, The non-linear redshift-space power spectrum: Omega from redshift 
surveys, MNRAS 279 (Mar., 1996) L1-L5, [astro-ph/9510049]. 

[61] E. Pajer and M. Zaldarriaga, On the renormalization of the effective field theory of large scale 
structures, JCAP 8 (Aug., 2013) 37, [arXiv: 1301.7182]. 

[62] A. Manzotti, M. Peloso, M. Pietroni, M. Viel, and F. Villaescusa-Navarro, A coarse grained 
perturbation theory for the Large Scale Structure, with cosmology and time independence in the 
UV, JCAP 9 (Sept., 2014) 47, [arXiv: 1407.1342]. 

[63] M. Pietroni, G. Mangano, N. Saviano, and M. Viel, Coarse-grained cosmological perturbation 
theory, JCAP 1 (Jan., 2012) 19, [arXiv: 1108.5203]. 

[64] M. Mirbabayi, F. Schmidt, and M. Zaldarriaga, Biased Tracers and Time Evolution, 
ArXiv:1412.5169 (Dec., 2014) [arXiv: 1412.5169]. 

[65] Z. Vlah, M. White, and A. Aviles, A Lagrangian effective field theory, ArXiv e-prints (June, 
2015) [arXiv: 1506. 0526]. 

[66] A. N. Taylor and A. J. S. Hamilton, Non-linear cosmological power spectra in real and redshift 
space, MNRAS 282 (Oct., 1996) 767-778, [astro-ph/9604020]. 

[67] M. Crocce and R. Scoccimarro, Nonlinear evolution of baryon acoustic oscillations, PRD 77 
(Jan., 2008) 023533, [arXiv:0704. 2783]. 

[68] T. Yano and N. Gouda, Evolution of the Power Spectrum and Self-Similarity in the Expanding 
One-dimensional Universe, ApJS 118 (Oct., 1998) 267-274, [astro-ph/9806026]. 

[69] A. E. Schulz, W. Dehnen, G. Jungman, and S. Tremaine, Gravitational collapse in one 
dimension, MNRAS ‘^31 (May, 2013) 49-62, [arXiv : 1206.0299]. 

[70] R. W. Hockney and J. W. Eastwood, Computer Simulation Using Particles. 1981. 

[71] G. Efstathiou, M. Davis, S. D. M. White, and C. S. Frenk, Numerical technigues for large 
cosmological N-body simulations, ApJS 57 (Feb., 1985) 241-260. 

[72] A. L. Melott, Massive neutrinos in large-scale gravitational clustering, ApJ 264 (Jan., 1983) 
59-78. 

[73] B. N. Miller and J.-L. Rouet, Cosmology in one dimension: fractal geometry, power spectra and 
correlation. Journal of Statistical Mechanics: Theory and Experiment 12 (Dec., 2010) 28, 
[arXiv: 1004. 0227]. 

[74] J. Carlson, M. White, and N. Padmanabhan, Critical look at cosmological perturbation theory 
techniques, PRD 80 (Aug., 2009) 043531, [arXiv:0905. 0479]. 

[75] P. J. E. Peebles, The large-scale structure of the universe. 1980. 


- 43 - 


